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In this chapter we review our knowledge of our galaxy’s cometary population outside 
our Oort Cloud - exocomets and Interstellar Objects (ISOs). We start with a brief overview 
of planetary system formation, viewed as a general process around stars. We then take a 
more detailed look at the creation and structure of exocometary belts, as revealed by the 
unprecedented combination of theoretical and observational advances in recent years. The 
existence and characteristics of individual exocomets orbiting other stars is summarized, before 
looking at the mechanisms by which they may be ejected into interstellar space. The discovery 
of the first two ISOs is then described, along with the surprising differences in their observed 
characteristics. We end by looking ahead to what advances may take place in the next decade. 


1. INTRODUCTION 


At the time of publication of Comets IT 
2004), the presence of cometary bodies around other stars 


had been well established. Since that time there has been a 
significant increase in theoretical investigations of cometary 
formation in protoplanetary disks (see §2.3]and Simon etal., 
this volume). Resolved imaging and spectroscopy of cir- 
cumstellar disks has become almost commonplace, giving 
exquisite insights into the birthplaces of comets (see §2). 
Increasingly common detections of gas in extrasolar Kuiper 
belts has opened a new approach to study exocometary gas 
and composition, combined with a new era of dust detection 
at IR and sub-mm wavelengths. At the same time, obser- 
vational studies have resulted in the discovery of new exo- 
comet host systems, the detection of dust continuum transits 
for the first time, and further understanding of extant sys- 
tems like 8 Pictoris (hereafter 6 Pic; see Fig. [I]and §3.2). 
White dwarf atmospheric modelling has revealed a variety 
of compositions of extrasolar small bodies, including po- 
tentially volatile-rich exocomets. 

Interstellar Objects (ISOs) are planetesimals unbound to 
stars, and generally thought to originate by ejection from 
their planetary systems by various physical mechanisms. 
These could be both asteroids that are generally consid- 
ered to be inert rocky bodies, or comets that contain a large 
amount of ices which can lead to outgassing and mass loss, 
the most common method of telling the difference between 
them (however, see the chapter by Jewitt & Hsieh in this 


Fig. 1.— Artist’s impression of the 8 Pic system, illustrat- 
ing the extensive families of exocomets observed within this 
system (credit ESO/L. Calçada). 


volume on the blurred boundaries between these objects). 
Assuming the same physics and dynamics occurs in other 
planetary systems as in our own, then the observable ISO 
population would primarily consist of cometary bodies of 
similar sizes to those orbiting our Sun, although see the dis- 
cussion in §6.4] 

In contrast to exocomets, the existence of ISOs was 
merely hypothesised in Comets IT, although their discov- 
ery was widely anticipated. Detection of ISOs passing 
through the Solar system would in principle allow remote 
(and eventually in-situ) sampling of bodies from other plan- 


etary systems. The surety that such a discovery would even- 
tually occur grew with increasingly refined models of plane- 
tary system evolution, and a growing understanding of how 
cometary bodies are lost to interstellar space and potential 
evolutionary processes, outlined in 

Increasing observational capabilities of wide-field sur- 
veys led to the first discovery of an ISO, II/‘Oumuamua, 
in 2017 (described in §6). That said, the second ISO 
discovery, 21/Borisov, was by an individual effort with a 
smaller telescope (see $6.2). 2I/Borisov was clearly an 
active comet from first detection, matching the expecta- 
tion that most ISOs would be exocometary in nature. But 
11/‘Oumuamua was very different, with initial observations 
failing to detect any sign of activity, even though it was ob- 
served at a heliocentric distance of 1.2—2.8 au. Hence, it is 
already clear that ISO properties span a wider range than 
previously thought, to the extent that it is difficult to predict 
what will be found in the next decade. 

In this chapter we will discuss the formation of and the 
observational evidence for exocomets, their ejection mech- 
anisms into space to become ISOs and the effects that this 
has on their physical properties, including processing in 
the interstellar medium. We then describe the discovery 
and characterization of the first two ISOs observed pass- 
ing through the solar system and the next decade of ISO 
science. 


2. PLANETARY SYSTEM EVOLUTION 


2.1. 


ISOs travel through the interstellar medium (ISM), 
which is the primary galactic repository from which stars, 
including their surrounding planetary system, form. Since 
the advent of 11/‘ʻOumuamua, it has become clear that the 
ISM contains, apart from gas and dust, a third component — 
interstellar objects. These are much smaller than the free- 
floating planets that have been detected via microlensing 
(Mróz et al.|2017) and directly imaged 
2022), but are much more numerous by many orders or 
magnitude. Thus ISOs are present during the entire process 
described in the following text. 

The star and planet formation process starts in the ISM 
when molecular clouds form through turbulent compression 
and global instabilities. Stars begin to form when denser 
parts of the clouds become unstable and start to collapse 
gravitationally. Once a molecular cloud begins to collapse, 
its central density increases considerably, eventually lead- 
ing to star formation, with often an entire star cluster emerg- 
ing (Bate et al.\2003). A large fraction of the formed stars 
are not single but binaries. The binary fraction is 30-50% 


for Solar-type stars in the local neighbourhood (Raghavan! 
et al.\2010) and up to 70% in young clusters hea 
van et al. 2010} Duchéne and Rraus|2013). This fact might 
be important because some of the suggested ISO ejection 
mechanisms rely on binaries being present (see 94.1). In 


this phase, a nascent individual protostar grows in mass via 
accretion from the infalling envelope until the available gas 


Star and primordial disk formation 


reservoir is exhausted or stellar feedback effects become 
important and remove the parental cocoon (for details, see 
and Bergin ef a. in this volume. 

As a natural consequence of angular momentum conser- 
vation, a disk develops during the formation of the protostar. 
Any initial nudge that imparts some core rotation means 
that material collapsing from its outer regions (with higher 
angular momentum) is channelled onto a disk, rather than 
the protostar itself (Zerebey et al.|1984). These young stel- 
lar objects (YSOs) consist of several components: the ac- 
creting protostar, the cloud (core) it is forming from, a disk 
surrounding it, and outflows of material that fail to be ac- 
creted either to the star or the disk. Observationally, YSOs 
are classified based on their continuum dust spectrum, thus 
the developmental stage of its disk. Class 1 objects are still 
deeply embedded within the cloud and mainly emit at mm 
and sub-mm wavelengths, leading to a characteristic posi- 
tive IR spectral index. Class 2 objects reveal sources with 
NIR and MIR excess, and Class 3 objects show only slight 
NIR excess. This classification scheme mainly reflects the 
development of the disks surrounding the stars from disk 
formation up to disk dispersal after several Myr. 


2.2. Disk properties and their development 


Flattened disks of cold dust and gas rotate around al- 
most all low-mass stars shortly after their birth 
jand Cieza\2011). These are the sites where planetesimals 
form. ISOs are generally believed to be ejected planetesi- 
mals. Therefore, we review the planetesimal formation pro- 
cess in order to understand better the properties of ISOs. 

Observations mainly at IR to millimeter wavelengths de- 
termine the frequency of disks and their mass, size, struc- 
ture, and composition as a function of age. The ages of indi- 
vidual young stars can only be determined with a relatively 
large error margin. Therefore, the temporal development of 
disks is usually determined in young stellar clusters where 
many relatively coeval stars are located in close proximity. 
The mean age of the cluster stars is taken to test the tempo- 
ral evolution of the disks in this environment. 

Disk masses seem to decrease rapidly with system age 
(Ansdell et al.\|2016). In the extremely early phases, the 
disk mass can be comparable to the mass of the proto- 
star, but after just 1 Myr, it is often less than 0.1% of the 
star’s mass. Some material is lost through outflows; some 
matter accretes onto the star or the star’s radiation photo- 
evaporates it. Significantly in this context, some material 
also condenses into centimeter- or larger-sized bodies, in- 
cluding planetesimals. It also seems that disk dust masses 
correlate with stellar mass (Raisk «x M2-°), but this corre- 
lation steepens with time, with a more substantial drop for 
low-mass stars Ansdell et al.{2016). 
Observations show that a surprisingly large fraction of disk 
dust masses appears low compared to the solid content of 


observed exoplanets (Najita and Kenyon||2014). This has 


been interpreted as a sign of early planet formation 
nara et al.\2018 Tychoniec et al./2020). However, recent 


work suggests that detection biases might play a role here 
and that disks contain similar amounts of solids as found in 
exoplanetary systems (Mulders et al.[2021). 

Observations suggest a shallower density profile in the 
inner disk (© œ r-t) and steeper at larger distances 
(ox r~?) (Andrews}|2020). However, the classical as- 
sumption of a smooth gas disk is not always appropriate. 
High-resolution observations with ALMA provide evidence 
for rings, gaps and spiral structures in at least some of the 
most massive disks. Pronounced gap structures are often 
interpreted as being carved by protoplanets, implying that 
planetesimal formation starts early in these environments. 
The current sample of very high-resolution measurements 
in the mm continuum or scattered light is biased in favor of 
larger, brighter disks that preferentially orbit more massive 
host stars Garufi et al.{2018). There- 
fore, any conclusions about the prevalence of substructure 
have to be carefully assessed. 

The typical disk sizes in dust measured by mm contin- 
uum observations lie in the range of 10 — 500 au. Gas disk 
sizes measured from the CO line emission are often larger 
than the dust disk size of the same disk. The disk dust radius 
appears to be correlated with the disk dust mass 
Sanchis et al. 2021). 

The disk lifetime is not only of prime importance for 
planet formation theory but equally for determining the 
timescale on which planetesimals form. The mean disk life- 
time is mostly determined from the observed disk frequency 
of clusters of different mean ages. It is found that the disk 
fraction decreases rapidly with cluster age with <10% of 
cluster stars retaining their disks for longer than 2-6 Myr 
in solar neighborhood clusters. However, recently several 
disks have been found with ages > 10 Myr, some as old 
as 40 Myr that still have the potential of forming planets 
(Flaherty et al./2019). Besides, the disk lifetime seems to 
be much shorter in dense clusters and longer in co-moving 
groups. The method of using disk fraction in clusters for 
determining disk lifetimes has been demonstrated to suf- 


fer from several biases (Pfalzner et al.|\2022). Looking 


at the transition phase from Class II to Class III objects, 
found a higher typical disk lifetime of 
x 8 Myr. A more complex picture is emerging, where disk 
lifetimes can range from <1 Myr to a few tens of Myr, and 
the mean age of the field star population lies in the range 
1-10 Myr, but this is still not definitely determined. 


2.3. Dust and planetesimal growth 


For planets that reside close to their star, dust accretion 
is the standard formation scenario. Here planet formation 
takes place over a variety of size scales (see Fig. B). Ini- 
tially, the mm-sized dust particles largely follow the move- 
ment of the gas. However, Brownian motion means that 
collisions between the dust particles take place neverthe- 
less. The low relative velocity of the dust particles means 
that they stick upon contact so that first larger porous frac- 


tal aggregates form. These compact in subsequent colli- 
sions and settle toward the midplane (Dubrulle et al.J1995), 
where the growth sequence continues. Once particles reach 
a size where they start to decouple from the gas aerodynam- 
ically, these solids start to migrate radially inwards. Further 
collisions lead to the formation of boulders, and eventually, 
planetesimals emerge. The largest planetesimals accumu- 
late nearby smaller ones, grow further, and become terres- 
trial planets or the cores of gas giants. However, it is the 
planetesimals that are of particular interest in this chapter. 
Today’s comets and interstellar objects are thought to be 
largely unchanged ancient planetesimals (for details, see, 


for example |Blum and Wurm||2008; Morbidelli and Ray- 
2016 2018} and Simon et al. of this book). 


This model has two potential problems: growth barriers 
during the accretion process and a relatively long formation 
timescale. The various growth barriers are the bouncing, 
fragmentation and drift barriers. The times expected for the 
different growth phases are 10? —10* yr to form mm- to cm- 
sized pebbles, 104 — 10° yr until the planetesimal stage is 
reached, 10°—10" yr to form terrestrial-type planets, and an 
additional 10° yr for the gas giants to accumulate their gas 
(Pollack et al.1996). The cumulative time seems at odds 
with observations of the disk frequency in young clusters, 
which indicate the median protoplanetary disk lifetime to be 
merely 1-3 Myr for both dust and gas (Haisch et al.)2001). 

However, over the last decade, mechanisms have been 
devised that might solve the timescale problem while at the 
same time overcoming the growth barriers (Armitage|2018). 
These mechanisms are the streaming instability (Youdin and 
Goodman\2005) and pebble accretion (Lambrechts and Jo- 
hhansen]2012). 

In the outer parts of the disks, planets also might form 
via a second process. In the initial phases, these disks are 
still very massive and can become gravitationally unstable. 
Simulations show that such disk instabilities lead to spiral 
arm formation, which can fragment directly to form large 
protoplanets (Boss||1997). However, the material near the 
star stays too hot to go unstable, so this process is expected 
to generate planets typically only at large distances (> 100 
au) from the star. Planetesimals would be produced as a by- 
product of planet formation. Details of the growth process 
including relevant references are described in the chapter on 
planetesimal formation by Simon et al. in this volume. 


2.4. Outer areas of young disks 


In we saw that the gas and dust disk differ quite of- 
ten in size. The difference in size is a result of the grain 
growth in the disk. Initially, the pressure gradient of a 
disk exerts an additional force that causes gas to orbit at a 
slightly subkeplerian speed. As long as the dust is coupled 
to the disk, it just follows the gas movement. However, as 
dust grains grow to about mm sizes, the orbiting grains ex- 
perience a frictional force, drifting inwards. At the same 
time, the disks’ viscosity means that the gas spreads out to 
conserve angular momentum and enable gas close in to ac- 
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Fig. 2.— Size scales during dust growth. Images from: (a) Chondritic interplanetary dust particle from Vessberger et al.| 
(2001) (CC license Attribution 2.5); (b) dust aggregate simulation from [Simon et al.) (2022); (c) 4.5 Gy old Allende 
meteorite with chondrules (CC Generic License 2.0); (d) images of boulders seen on by Rosetta comet 67P on 16 May 
2016 from a distance of 8.9 km (ESA/Rosetta/NAVCAM - CC BY-SA IGO 3.0) (e) nucleus of 67P seen by Rosetta on 7 
July 2015 from a distance of 154 km with a resolution of 13.1 m/pixel (ES A/Rosetta/NAVCAM - CC BY-SA IGO 3.0). 


crete onto the star. Observations at (sub-)mm wavelengths 
typically trace the large dust grains; thus, disks appear more 
extended in gas than dust. The combined effects of growth 
and vertical and radial migration of dust particles mean that 
the disks’ mean particle size should decrease with distance 
to the midplane and the central star. Alternatively, in terms 
of timescales, dust growth should proceed slower off the 
midplane and in the outer disk areas (Dullemond and Do- 
where do the planetesimals that are ejected and become ISO 
primarily form? 

In the outer disk areas, at particular distances from the 
central protostar, the temperature is so low that volatile 
compounds such as water, ammonia, methane, carbon 
dioxide, and carbon monoxide can condense into solid ice 
grains. These distances are referblack to as the snowlines. 
These snowlines are essential in the context of dust growth. 
Ices can change the effective particle strengths and thereby 
affect collision outcomes. argue that the 
critical velocity for fragmentation increases at the (CO or) 
CO» and NH3 snowlines and decreases at the H2O snow- 
line. 

However, direct imaging has observed exoplanets at 
large distances (> several tens of au) from their parent stars. 
Either these planets have formed there or migrated outwards 
to these orbits. Among these directly imaged planet hosts, 
there are several that are thought to be relatively young 
(< 20 Myr; for example, PDS 70 (Mesa et al.[2019)), which 
illustrates that even in the outer disks, planetesimals must 
be able to grow relatively fast. If protoplanets can grow on 
such a timescale, even in the outer disk, planetesimals and 
the resulting ISOs should form on even shorter timescales. 


3. EXOCOMETS 


Exocomets are small bodies orbiting other stars that ex- 
hibit signs of activity, through the release of dust and/or 
gas (Strøm et al.{2020). As such, their definition is broader 
compablack to the solar system population, which is much 
more refined dynamically and compositionally. Evidence 


for exocomets dates back to 1980s, when exocomets were 
inferblack from variable absorption features in the spectrum 
of 6 Pic (Ferlet et al.[1987), and when the dust from an exo- 
cometary belt was first discoveblack around Vega 
et al.|1984). In general, exocometary release of gas and dust 
is detected in both the inner (few stellar radii to few au) and 
the outer (~ 10 to 100s au) regions of extrasolar planetary 
systems (as depicted in Fig. Bp. and we will distinguish be- 
tween these two populations in the following sections. 


3.1. Exocometary belts: extrasolar cometary reser- 
voirs 
3.1.1. Detection and basic observational properties 


Exocometary belts (also known as debris disks) are ex- 
trasolar analogs of our Kuiper belt, and are the reservoirs 
of icy exocomets in planetary systems. These belts are first 
detected through the dust produced by exocomets as they 
collide and grind down to release small, observable dust. 
Thermal emission as well as scattered starlight by this dust 
can be detected and used to study these belts. The major- 
ity of detectable belts lie at distances of tens of au from 
their host star, and are cold with typical temperatures of 
tens of K (Chen et al.|[2014). They are therefore first de- 
tected by mid- to far-IR surveys, most sensitive to the peak 
of the dust’s thermal emission, with space telescopes such 
as IRAS, Spitzer, WISE, and Herschel (e.g. {Su et al.[2006} 


Sibthorpe et al.|2018). These surveys show belts to be ubiq- 
uitous; they are detected around ~20% of nearby, several 


Gyr-old field stars, even though the surveys’ sensitivity at 
best enabled detection of belts a few times more massive 
than the Kuiper belt (Eiroa et al.|/2013). This makes this 
number very much a lower limit; studies of younger, less 
evolved belts at a few tens of Myr indeed show occurrences 
as high as ~75% (Pawellek et al.[2021). 

The typical observable (i.e. not the total) mass of belts, 
only considering solids of sizes up to ~cm in diameter, 
range between ~0.1-25% of an Earth mass, or about a tenth 


to 20 times the mass of the Moon (Holland et al.\\2017). 


This is but a small fraction of the belts’ total mass, which 
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Fig. 3.— Top-down view of a young, ~ 10 — 100 Myr-old planetary system after dissipation of the protoplanetary disk. 
The system has largely formed outer gas and ice giant planet(s) (at few to tens of au, like the Neptune analog with its 
orbit depicted as the dashed line), has an exocometary belt at tens of au producing cold gas and dust, and star-grazing 
exocomets in the inner region, producing gas and potentially warm exozodiacal dust as they approach and recede from the 
star. Exocomets can be probed in a variety of ways. In the outer regions (blue/cyan symbols), cold exocomets (extrasolar 
KBO analogs) can be probed in the belts they formed in, through cold dust and gas emission, and cold absorption against 
the star for edge-on systems. In the inner regions (black symbols), exocomets are probed through warm exozodiacal dust 
emission, as well as blue- and black-shifted absorption (see also Fig. |5) and asymmetric ’shark-fin’ transits for edge-on 


systems. 


is dominated by the largest, unobservable bodies. The to- 
tal mass of the Kuiper belt is estimated from observations 
of its largest bodies, and is about 6% of an Earth mass 
[Ruscio et al.|2020). The Kuiper belt’s dust mass is how- 
ever predicted to be much lower than that of exocometary 
belts, to the point that a true Kuiper belt analog would not 
be detectable around nearby stars (Vitense et al.[2012}. In 
summary, detectable exocometary belts have dust masses 
typically ~ 10 — 10° times higher than inferred for the 
present day Kuiper belt, where this broad range represents 
the spread of belt dust masses observed. 


3.1.2. Birth of exocometary belts 


Exocometary belts are born at tens to hundreds of au 
from the central star within protoplanetary disks, but they 
cannot be identified until after the protoplanetary disk is 


dispersed. During dispersal, primordial gas and dust are 
removed leaving behind second-generation dust (and gas), 
produced destructively in collisions of larger bodies. The 
dispersal of the protoplanetary disk takes place rapidly once 
the star is a few to ~10 Myr-old 


2017| and §2.2). The emergence of an exocometary belt 
from this dispersal is evidenced by the presence of an in- 


frared excess, but with a decrease of ~2 orders of magni- 
tude in dust mass compared to proto- 
planetary levels. 

In practice, the largest exocomets within belts must form 
during the protoplanetary phase to access the gas needed 
for the planetesimal formation process (see the chapter by 
Simon et al., this volume), but observationally validating 
this is difficult. An interesting proposal is that exocometary 
belts may originate in (one or more of) the large rings 


observed in structured protoplanetary disks (Michel et al. 
2021), which has been shown could potentially explain 


the population of observed exocometary belts (Najita et al.| 
2022). Another possibility is that the photoevaporative dis- 
persal of the gas-rich protoplanetary disk induces high dust- 


to-gas ratios (Throop and Bally}2005) that may trigger the 


streaming instability and produce massive belts of planetes- 
imals/(exo)comets (Carrera et al.|2017). However, such 
massive, large belts are not observed, and other models of 
dust evolution in a photoevaporating disk indicate that the 
grains efficiently drift inwards faster than they can pile-up 
to trigger planetesimal formation (Sellek et al.|2020). 

One difficulty in the protoplanetary disk - exocometary 
belt observational comparison is that most of the protoplan- 
etary disks in nearby star-forming regions surround low- 
mass stars, whereas the majority of detectable exocometary 
belts orbit around A-F stars, likely due to current lack of 
sensitivity to belts around later-type stars, particularly M 
dwarfs (Luppe et al.|/2020). In general, the presence of 
evolved Class III disks in young star forming regions, with 
masses consistent with young exocometary belts, suggests 
that exocomets can form early on, within the first ~2 Myr 
(Lovell et at {2021}. 

A key component to the transition from a protoplanetary 
disk to an exocometary belt is the presence of gas. This 
produces drag, slowing collisions and favoring dust growth 
in protoplanetary disks; on the other hand, gas dispersal en- 
ables the more collisionally destructive, dust-producing en- 
vironment of exocometary belts (Wyatt et al.[2015). Unfor- 
tunately, distinguishing primordial from second generation 
gas is even harder, due to the surprisingly low CO abun- 
dances observed in many protoplanetary disks 
[2020}, combined with the presence of gas now discovered 
in a growing number of exocometary belts (see 3.1.5), and 
the difficulty in measuring total gas masses in both evolu- 
tionary phases. 


3.1.3. Evolution of exocometary belts 


The physics of observable (massive) exocometary belts 
is driven by collisions, causing the grind-down of large ex- 
ocomets (analogs to the solar system’s observable Kuiper 
belt objects) into small, observable grains. This produces 
a collisional cascade, which rapidly reaches a quasi-steady 
state size distribution where the number of solid particles 
scales as n(a) x a1, where q = —3.5 for a cascade extend- 
ing to infinitely large and small bodies (Dohnanyi\|1969). 
More detailed cascade simulations show that this q value 
can vary between -3 and -4 
[Schlichting|2012). Observations of mm-cm wave spectral 
slopes implying values varying between ~ — (2.9 — 3.8) for 
observable grains in the mm-cm size regime (Norfolk et al.) 
(2021), although interpretation of measurements in terms of 
size distributions is complex (Léhne|[2020). These slopes 
imply that what extrasolar observations are sensitive to (the 
solids’ cross-sectional area) is dominated by small up to 
~cm-sized grains, but the total mass is dominated by the 


largest, unobservable bodies. This is the opposite to obser- 
vations of the Kuiper belt, where we are most sensitive to 
the largest bodies. 

At the bottom end of the cascade, grains are sensitive 
to radiation forces, and are typically removed by radiation 


pressure from the central star (Burns et al.|1979), or by stel- 


lar winds which are particularly important for young late- 
type stars 
2006). For a radiation pressure-dominated environment, the 
smallest, blow-out grain size in the cascade is the one where 
the radiation pressure force dominates over gravity, pushing 
smaller grains on unbound orbits and rapidly ejecting them 
from the system (Backman and Paresce|1993). 

This removal of small grains causes a net mass loss of 
material from the cascade as a function of time. In the 
simplest model of a collisional cascade, this would cause 
a belt’s mass to initially stay constant - until the largest 
bodies begin colliding - and then decrease as t7! 


2003) although models with more real- 


istic treatment of cascade processes predict somewhat dif- 
ferent time decays 
(2017). Nevertheless, a decrease in dust emission as a func- 
tion of time is unambiguously observed in surveys of de- 
bris disks of different ages, and can largely be explained by 
collisional evolution models (Sibthorpe et al.|2018). Addi- 
tional mass removal effects could come into play that may 
speed up this mass loss, as hinted at by recent observations 
(Pawellek et al.{2021). Notably, faster mass removal could 
be caused by gravitational interaction between planets and 
exocometary belts, potentially ejecting small bodies from 
the planetary system and contributing to the population of 
interstellar objects (see f. Such interaction could be akin 
to the late dynamical instability and/or to the outward mi- 
gration of Neptune inferred to have depleted the vast major- 


ity of material in our solar system’s Kuiper belt (Malhotra 
1993!|Gomes et al.|2005). 


3.1.4. Resolved dust imaging: structure and (exo)planet- 
belt interaction 


Resolved imaging is ultimately needed to study the lo- 
cation and structure of exocometary belts. Given the star’s 
luminosity, the observed temperature from the dust’s unre- 
solved spectrum can be used to calculate its location/radius 
(distance from the star), under the assumption that the dust 
behaves as a blackbody. In reality, small dust dominating 
the spectrum’s emission is an inefficient emitter, and thus 
such temperature-derived blackbody radii typically under- 
estimate the true belt locations by a factor of a few 
feral 2073). 

Belts are detected around stars from as far as nearby star- 
forming regions (~ 150 pc) to our nearest neighbors (a few 
pc away), with typical angular sizes of a few tenths to a few 
tens of arcseconds on-sky, so easily resolvable with the lat- 
est generation of high resolution facilities from optical to 
mm wavelengths. Given their low masses and large on-sky 
extents, the difficulty in imaging exocometary belts arises 
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Fig. 4.— Images of debris disks and exocometary belts. Fomalhaut credit: ALMA (ESO/NAOJ/NRAO), M. MacGregor; 
NASA/ESA Hubble, P. Kalas, B. Saxton (NRAO/AUI/NSF); other images are adapted from{Marino|(2022). 


often from their low surface brightness, requiring deep in- 
tegrations. Nevertheless, several famous belts are bright 
enough for structural characterization at high resolution. 

The first resolved exocometary belt was that around 8 
Pic (Smith and Terrile|1984). Almost four decades later we 
now have a considerable inventory of images for a few tens 
of belts (e.g. Fig. (4). These have been obtained across 
the wavelength spectrum, from space and ground-based 
optical/near-IR scattered light images (milliarcsec resolu- 
tion, e.g. 
et al.\2020), to space-based IR (few arcsec resolution, e.g. 
and mm-wave ground-based interfero- 
metric images (sub-arcsecond resolution, e.g.[Marino et al.| 
Sepulveda et al./2019). 


The most basic measurable quantity of a resolved belt is 
its location (radius). Exocometary belts resolved across the 
wavelength range show typical radii of a few tens to 100s 
of au, with a shallow trend indicating that more luminous 


stars tend to host larger belts (Matra et al.\2018a\|Marshall 
2021). While a combination of observational bias 


and radius-dependent collisional evolution can explain this 
positive trend, the low scatter suggests that belts (includ- 
ing our own Kuiper belt) may form at preferential locations 
in protoplanetary disks. If confirmed with a larger sample 
of belts, this would imply a potential connection to either 
temperature-dependent exocomet formation processes (e.g. 
the CO ice line), or to the timescale of planet formation 
processes preventing exocomets from growing further into 
planets at these distances. 

While the classic picture of an exocometary belts is that 
of a narrow ring, recent imaging at IR and mm wave- 
lengths shows that broad belts (with width/radius, AR/R, 
approaching ~ 1) are at least as common as narrow rings 
(AR/R « 1; see Fig. |4). Some of these belts may be born 
narrow but others, like the famous HR 8799 belt, may have 


evolved to host a scattered disk of high eccentricity particles 
akin to our Kuiper belt’s, through gravitational interaction 
with interior planets (Geiler et al.[2019). 

Beyond a simple radius and width, exocometary belts 
observed at sufficiently high resolution show structural fea- 
tures which have long been linked to the gravitational ac- 
tion of planets (or dwarf planets) interior, exterior, or within 
these belts (e.g.|Pearce et al.|2022). (Sub-)structures can be 
broadly separated into radial, vertical, and azimuthal fea- 
tures. 

In the radial direction, a belt’s inner/outer edge and its 
sharpness can be linked to dynamical truncation by plan- 
ets just interior or exterior to it (Chiang et al.|2009). This 
is particularly true for inner edges with surface density pro- 
files steeper than the r“ with a ~ 2—2.3 expected from col- 


lisional evolution of an undisturbed disk (Schiippler et al. 
2016; |Matrà et al.\2020). 


In the absence of truncating planets, the sharpness of belt 
edges may also be used to probe the eccentricity dispersion, 
and hence the level of dynamical excitation of exocomets 
within the belt [Rafikov||2023). Dynamical 
excitation of exocometary material within belts can also be 
probed by the belts’ observed vertical structure in edge-on 
systems [Daley et al.{2019). Like the 
Kuiper belt, this has the imprint of the system’s dynamical 
history; for example, the vertical structure of the 6 Pic belt 
supports the presence of a double population of exocomet 


inclinations (Matra et al.|2019b), akin to the hot and cold 
populations of Kuiper belt objects (KBOs) in the solar sys- 


tem (Brown|2001). 


While the radial distribution of some wide exocometary 
belts appears to be smooth 
et al.|2021), at least three systems so far (e.g. |Marino et al. 
2018b) show the presence of rings separated by a gap, rem- 
iniscent of the rings observed in younger protoplanetary 


disks (Andrews et al.|/2018). In the simplest scenario, a 


single stationary planet at the gap location (tens of au) can 
clear its chaotic zone and carve a gap of width and depth 
related to the planet’s mass and the age of the system 
rison and Malhotra\2015). However, these gaps may be too 


wide and require planets to migrate (Friebe et al.|2022), or 


different scenarios such as secular effects due to an eccen- 


tric planet interior to both rings (Pearce and Wyatt|2015). 
A few narrow belts, most prominently the famous Fo- 

malhaut and 

HD202628 rings 


are observed to be significantly eccentric (e ~ 0.12). Sec- 
ular perturbation by eccentric planets interior to the belts 
can produce and maintain eccentric belts (e.g. 
[1999)}, though the narrower than predicted width of these 
rings may point to alternative scenarios, such as belts be- 
ing born eccentric in the protoplanetary phase of evolution 
(Kennedy|2020). 

Azimuthal asymmetries in the brightness of exocometary 
belts are often detected in scattered light observations, par- 
ticularly with the Hubble Space Telescope (HST; 
der et al.\2014), though these may be simply produced by 
a combination of disk eccentricity, grain scattering phase 
functions, and viewing geometry (Lee and Chiang||2016), 
rather than true azimuthal asymmetries in their density dis- 
tribution. Similarly, eccentric disks observed in thermal 
emission may produce brightness asymmetries due to peri- 


center and/or apocenter glow (e.g.|MacGregor et al.|2022), 
which which are expected to vary with wavelength 


[Eynch and Lovell|2022). The strongest evi- 
dence for azimuthal asymmetry arises from the edge-on 8 
Pic belt, which displays a strong clump in dust emission at 
mid-IR wavelengths (Telesco et al.12005 2023) 
but not as prominently in scattered light or ALMA dust 


observations (Apai et al./2015 Matra et al.J2019b). Such 


clumps may be produced by resonant trapping by a migrat- 


ing planet at the inner edge of the belt (Wyatt|2003), akin to 


Neptune’s migration into the Kuiper belt in the young solar 
system (Malhotra\1993), but also by giant impacts between 
a a or aT Sty) 

With the clear dynamical link between exocometary belt 
structure and the presence of exoplanets interacting with 
them, it is natural to ask whether the very existence of a belt 
is linked to the presence (or not) of exoplanets in the same 
system. For example, if the brightest debris disks form from 
the most massive protoplanetary disks that formed exoplan- 
ets most efficiently, one may expect a positive correlation 
between the presence of a belt and planets 
[2007). On the other hand, the presence of outer giant plan- 
ets may efficiently eject and deplete exocometary belts (pro- 
ducing ISOs), while enabling terrestrial planet formation to 
take place interior to them without significant disruption; in 
this case, one would expect an anticorrelation with outer gi- 
ant planets, and a correlation with low-mass planets 
(2012). Evidence of trends continues to 


be debated in the observational literature: a correlation be- 


tween low-mass planets and belts (through their infrared ex- 
cess) had initially been found 
et al.|2014), but a study with an expanded sample found no 
correlations with either low- or high-mass planets 
Martin et al.\2015). Similarly, no correlation was found be- 
tween the presence of radial velocity planets and the prop- 
erties of exocometary belts (Yelverton et al.|/2020). It is 
notable however that a significant number of long-period, 
young directly imaged planets are found in systems with 
bright debris disks, with some evidence of a correlation 
(Meshkat et al POTN). 

In summary, the structure of exocometary belts at tens 
of au probes the orbits of the icy bodies within it, and, as 
is the case in our solar system’s Kuiper belt, can reveal dy- 
namical interactions with mostly unseen outer planets. As 
proven by dynamical simulations successfully reproducing 
these structures, these gravitational interactions can remove 
exocomets in two ways. They can either eject exocomets 
from outer belts, potentially producing Oort clouds analogs, 
or directly turning them into ISOs (§4.1.3), but also to scat- 
ter them inwards, which could eventually produced the tran- 
siting exocomets observed close to the central stars ($3.2) 
and deliver volatile-rich material to planets in the habitable 
zones of their host stars. 


3.1.5. Evidence for gas within outer exocometary belts 


Since their discovery almost 4 decades ago, exocometary 
belts have been considered to be gas free, to the point 
that this was long considered a defining difference be- 
tween younger protoplanetary and more evolved debris 


disks (Zuckerman et al.}1995). However, evidence for cir- 


cumstellar gas in these systems was discovered as early as 
1975 around ( Pic (Slettebak||1975), in the form of nar- 
row Ca II H and K absorption at the stellar velocity (i.e. 
not significantly red- or blue-shifted, as is instead treated 
in §3.2) which was later spatially resolved to extend out to 


100s of au from the star, where the dusty belt lies(Brandeker| 
2004). For a long time, evidence for cold gas in these 


systems remained limited to the 6 Pic and 49 Ceti belts. 
Their near edge-on geometry allowed detection of gas in 
both emission and absorption against the central star (Zuck- 
2014). 


The evidence for exocometary gas has rapidly strength- 
ened in the past decade. The sensitivity advance to cold 
gas emission brought by Herschel and especially ALMA 
has led to detections of atomic (C I, O I and/or C II) and/or 
molecular (CO) emission within over 20 exocometary belts, 
largely around ~10-40 Myr old stars (Moor et al.|/2017), 
but also as old as 440 Myr a ET, 
Early ALMA CO gas detections highlighted a dichotomy 
between belts with very high CO masses, approaching those 
of protoplanetary disks (107? — 107! Mg, e.g. 
2013), and more tenuous CO gas disks (1077 — 107® Ma, 
e.g.|Matra et al.201Ta). 


The origin of the gas remains an open question for high- 


mass CO systems. In a primordial scenario, the gas is a 
remnant of the protoplanetary disk (contrary to the dust, 
which is second-generation, leading to the term hybrid 
disk); Hə dominates the gas mass, allowing CO to remain 
shielded from stellar and interstellar photodissociation over 
the lifetime of the system (Kóspál et al./2013). In a sec- 
ondary scenario, CO gas is continuously produced by ex- 
ocomets within the belt, just like the dust (Zuckerman and! 
Marino et al 2018). How- 
ever, the production rate is unlikely to be sufficient to main- 
tain the high CO masses observed, so a shielding agent 
other than Hə (which is only produced in trace amounts 
in solar system comets) is needed. Shielding CO pho- 
todissociation (producing C and O) by C photoionization 
is a promising scenario to produce what would be second- 
generation, shielded high-mass CO disks (Kral et al.[2019). 

For low-mass CO belts, like Fomalhaut, HD181327 and 
B Pic, the secondary scenario is heavily favored, as the CO 
photodissociation timescale, even when accounting for the 
potential of unseen Hə with primordial CO/H2 ratios, is 
much shorter than the age of the system (e.g. 
(2016). This simple argument implies that CO must have 
been recently (and likely continuously) produced, and must 
be of secondary, exocometary origin. 


3.1.6. Gas release processes, composition, and evolution 


The gas release process of exocomets within belts from 
their host stars is likely different from solar system comets. 
In the very cold (~tens of K) environments at tens of au, 
the heating and sublimation-driven process producing gas 
in single solar system comets and ISOs (as they approach 
the Sun) is unlikely to be efficient; while activity sug- 
gesting CO outgassing takes place out to Kuiper Belt dis- 
tances in at least some long period comets (e.g. C/2017 


K2 (Pan-STARRS) 2017a), no activity is de- 
tected from KBOs on stable orbits (Jewitt et al. (2008) in- 
cluding Arrokoth (486958) 2014 MU69, |Lisse et al.|2021). 


Care should be taken in associating locations with temper- 
atures in the Kuiper belt versus exocometary belts. As ex- 
ocometary belts are most commonly observed around more 
luminous, A-F type stars, they will have warmer tempera- 
tures at the same radial location. This is however mitigated 
by the temperature dependence with stellar luminosity be- 
ing shallow for blackbody-like grains heated by starlight 
(T x L°-?5), combined with the fact that belts tend to be 
located at larger radii around more luminous stars (§3.1.4). 

One possibility is that we are observing a vigorous sub- 
limation period for exocomets immediately following the 
dispersal of the protoplanetary disk 
Kral et al.\2021). In this period, hypervolatiles like CO, N2 
and CH, sublimate at high rates due to the heating by newly 
visible starlight, previously shielded by large amounts of 
protoplanetary dust. In this scenario, Arrokoth-sized ob- 
jects (~10s of km diameter) at ~45 au from the central star 
take 10-100 Myr for the heat to reach their interior, which 
could explain the higher detection rate around young exo- 


cometary belts in that age range. However, in exocometary 
belts a large range of sizes would have to be considered, 
and it remains unclear whether the production rates could 
be sufficiently high to explain the observed CO masses. 

Other physical mechanisms have been proposed to re- 
lease gas in the collisionally active environment of young 
exocometary belts. While typical collision velocities are 
unlikely to result in sufficient heating, collisional vaporiza- 
tion of very small grains accelerated by radiation pressure 
could lead to ice release into the gas 
(2007). Additionally, UV photodesorption could very effec- 
tively remove volatiles such as H2O from the surface of icy 
grains/objects (Grigorieva et al.)2007), given these are con- 
tinuously produced/resurfaced through the collisional cas- 
cade. 

Once released, in the absence of shielding, molecular gas 
is expected to be rapidly destroyed by the stellar and inter- 
stellar UV fields. CO (with N2) is the most photo-resistant 
molecule, which may explain why it is the only one de- 
tected so far (Matra et al.[2018b). The CO ice mass fraction 
of exocomets can then be estimated from the ALMA mass 
measurements and steady state production/photodestruction 
arguments, indicating an exocometary CO content within 
an order of magnitude of solar system comets 
- showing that accessing exocometary 
compositions similar to solar system comets is now possi- 
ble. Unfortunately, the short photodissociation timescales 
and/or weak transitions of other molecules make other 
molecular detections currently challenging 
BOI8b|[KTusmeyer et al.2021). 

On the other hand, atomic gas is abundant and detected 
both in emission and absorption for edge-on belts. The de- 


tection of atomic N (Wilson et al.\/2019) and S, as well as 
C and O (see [Roberge et al.|2006\ and references therein) 


in the 8 Pic belt suggest that molecules other than CO are 
being released, arguing against the gas release being lim- 
ited to hypervolatile molecules. Atomic gas accumulating 
over time at the belt location is expected to produce a disk 
that viscously spreads over time, at a rate determined by the 
poorly-constrained gas viscosity, and eventually form an ac- 
cretion disk (Kral et al.\2016). This model can largely ex- 
plain the population of CO masses observed 
2020), including high CO mass disks through C shield- 
ing. However, resolved observations of C I gas by ALMA 
show inner holes inconsistent with accretion disk profiles, 


suggesting a more complex evolution (Cataldi et al.\2018 
2020). 


3.1.7. Detectability of exo-Oort clouds 


Extrasolar Oort cloud analogs may also exist, and ex- 
ocomets in those stars would be heated by absorption of 
radiation from their central star, nearby stars and the Cos- 
mic Microwave Background (CMB). The challenge for de- 
tecting extrasolar Oort clouds is their low surface density, 
low temperature and large sky areas, which is why long- 
wavelength surveys around nearby stars are best suited for 


detection. An additional factor is whether dust might be 
produced in Oort clouds, and if so how long it might sur- 


vive. (1991) used IRAS to search for excess 


far-IR emission around 17 nearby stars but did not achieve 
any detections. More recently, used 
Planck survey data to average the sub-mm emission around 
selected stars within 300 pc, and also performed a directed 
search around Fomalhaut and Vega, but again did not detect 
any signal. Intriguingly, when stacking 43 hot stars within 
40-80 pc, they detected a statistically significant excess flux 
at radial distances of 10* — 10° au. Searching around early- 
type stars is a key factor, as this will increase stellar heat- 
ing of the exo-Oort body surfaces and aid detection. But 
they cautioned that this could be a false positive due to the 
fine structure of galactic dust on scales of the Planck beam 
width. Alternatively it could be a real signal caused by weak 
nebular emission or dust grains ejected from debris disks 
near the stars. Future CMB missions such as LiteBIRD 
(Hazumi et al.\2020) may approach the sensitivity required 
for unambiguous detection. However, the 0.5° mission res- 
olution would restrict it to searches around the nearest stars 
within 10 pc, for which 104 au spans ~ 0.3°. Unresolved 
emission from more distant stars would be ambiguous as it 
could also be caused by thermal emission from dust closer 


to the star, see §3.1.4| 


3.2. Exocomets in inner planetary systems 
3.2.1. The P Pictoris system 


The star 6 Pic is a nearby (d = 19.3 pc) A6V star that 
was the first to have a circumstellar debris disk optically 
imaged (Smith and Terrile|1984). Soon after, transient red- 
shifted absorption features (see Fig. [5) in the UV spectrum 
were detected in the wings of the stellar Al II, Mg II and 
Fe II lines (Lagrange-Henri et al.J1988), and in numerous 


high-resolution optical spectra that include the strong Ca II 
H&K lines (Ferlet et al.|/1987). These were interpreted 
as originating from planetesimals undergoing sublimation 
when within a few stellar radii from 8 Pic, with most fea- 
tures being redshifted from motion towards the star. They 
were subsequently termed ‘Falling Evaporating Bodies’, or 
FEBs. In modern parlance these are exocomets, but care 
must be taken. The observation of metal lines may imply 
sublimation of refractory elements, and hence not neces- 
sarily a volatile-rich body as is common in our solar sys- 
tem. Subsequent studies of these exocomets have tended to 
use the strong Ca II H & K lines in high-resolution optical 
spectra, although Na I D-line absorption is also seen. Un- 
fortunately only metal lines in these exocomets have been 
measured and the lack of detections of the volatile species 
seen in solar system comets has prevented a direct compari- 
son. While studies of the gas in the outer circumstellar disk 
show it is likely to be released by exocomets in those re- 
gions as discussed in there is as yet no direct link to 
the exocomets seen close to the star. 

Akin to solar system comets, these exocomets are clearly 
undergoing significant mass loss, which argues for a source 
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Fig. 5.— Spectra of the core of the Ca II K stellar line 
in 8 Pic at three separate epochs over 4 days. The stable 
central narrow absorption is due to gas in the circumstellar 
disk, marking the stellar radial velocity. Transient Ca II 
absorption features caused by individual exocomets can be 
seen both blue-shifted and red-shifted relative to the star 
(see and Fig. B). Data from the ESO 3.6-m+HARPS 
public archive (see acknowledgments). 


region in the system. |Beust and Morbidelli| (2000) and 
Thébault and Beust| (2001) had previously shown that a 


Jupiter-sized planet within 20 au could sustain such a pop- 
ulation by strong orbital evolution via the 4:1 and/or 3:1 
mean-motion resonances. This would place the source re- 
gion at 4-5 au from £ Pic, within the imaged debris disk. 
The discovery by direct imaging of the Jovian exoplanet 
B Pic b by [Lagrange et al.| (2009) orbiting ata = 9.9 au 
gave credence to this hypothesis (see above). The 
more recent discovery of the closer massive planet 8 Pic c at 
a = 2.7 au will result in more com- 
plex dynamical evolution, but clearly enhances the possi- 
bility of scattering exocomets existing in the observed disk 
into star-approaching orbits. 

Although exocomets typically transit the star over a cou- 
ple of hours, was able to measure radial 
acceleration by individual bodies, and constrain their or- 
bital parameters via MCMC methods. A larger analysis 
was made by of 252 individual exo- 
comets seen in ~ 6,000 spectra obtained in 2003-2011. 
The velocities and absorption strengths display a bimodal 
distribution, implying at least two dynamical populations. 
Exocomets creating shallow absorption lines tend to have 
small total absorption and higher velocities, with perias- 
tron distances ~ 0.08 au (9 stellar radii) over a wide range 
of longitudes. The less numerous exocomets that produce 
shallow absorption also exhibit smaller radial velocities due 
to larger periastron distances of ~ 0.16 au (18 stellar radii). 
A small range of longitudes of periastron imply these less 
active objects all share similar orbits, and may result from 


the breakup of a single progenitor similar to the sungrazing 
Kreutz family (Vones et al.|2018]and references therein). 
Exocomets orbiting 8 Pic have now also been detected 
via broadband optical light curves using the Transiting Exo- 
planet Survey Satellite (TESS) mission (Zieba et al.[2019). 
The transit signals from individual exocomets show a ‘shark 
fin’ appearance caused by the asymmetric morphology of 


dust comae and tails, as predicted by|Lecavelier des Etangs 
(1999). Analysing all TESS light curves at that time, 


identified 30 individ- 
ual transits. Assuming that each exocomet on average has 
the same fractional surface area sublimating, the differential 
size distribution has a power-law exponent of —3.6 + 0.8, 
similar to that seen in solar system comets. 


3.2.2. Other systems exhibiting exocomet signatures 


As explored in Section[2| the decreasing mass of circum- 
stellar material as a young star begins its main-sequence 
lifetime implies that exocomets may also be more abun- 
dant at early times. Indeed, transient spectroscopic absorp- 
tion features similar to Pic are seen around many Herbig 
Ae/Be stars (Grady et al.|1996). With £ Pic as the archetype 
main-sequence host star, exocomets have been confirmed 
via similar spectroscopic or photometric transient features 
around 3 other stars: 49 Cet, HD 172555 and KIC 3542116 
(Strøm et al.\2020|and references therein). In addition there 
are likely exocomet systems such as c Aql (HD 183324), 
which exhibits Ca II variability on timescales of minutes 
[elesias er al](2018). In 
total, spectroscopic and photometric signatures indicating 
probable exocomet activity have been reported for another 
30 or more stars (Strøm et al.[2020). 

The common factor in all of these stars is that they are 
early-type A-F stars with bright apparent magnitudes. The 
brightness allows high signal-to-noise photometry from 
space-based missions such as Kepler and TESS to reveal 
weak transit signals from stars with less stellar activity than 
later K/M-type stars, while B/A stars also exhibit relatively 
clean optical spectra with fewer stellar absorption lines, 
facilitating spectroscopic detection. In terms of total num- 
bers, a directed spectroscopic survey of 117 B-G stars by 
found ~ 15% exhibited transient 
Ca II and Na I absorption. In summary, there is strong 
evidence that exocomets are commonly associated with 
early-type B/A stars. The evidence for later-type stars is 
less extensive, given that a search of 200,000 Kepler light 


curves by (2019) only found 5 potential pho- 


tometric transits occurring for 3 stars. 


3.2.3.. Warm exozodiacal dust 


Aside from the more distant cometary belts described in 
approximately 20% of A-K stars exhibit thermal sig- 
natures of warm (T ~ 300 K) dust within a few au (see|Kral| 
[et al.|2018). The combination of the Poynting-Robertson 
effect on large dust grains and radiation pressure on small 
dust grains implies very short lifetimes without continuous 
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resupply. Potential sources of this dust are either asteroidal 
collisions or dust ejection from comets. This may espe- 
cially be the case for systems with colder outer dust disks, 
where the source of the dust is expected to be ice-rich bod- 
ies. Given that the generation of the solar system’s zodiacal 


cloud is dominated by short-period comets (Rigley and Wy- 
|att|2022), it is plausible that this may be also the case for 


some exosystems (Marboeuf et al.[2016). A survey of exo- 
zodiacal dust around 38 stars showed a significant correla- 
tion with the existence of cold dust disks but no correlations 
with stellar age, as would be expected by analogy with the 
solar system (Ertel et al.|2020). 

An exemplar star where this appears to be the case is the 
main-sequence F2 star 7 Corvi. This star possesses a cold 
outer belt centered at ~ 165 au, plus a warm inner dust ring 


whose inner boundary is as close as 2-3 au (Duchéne et al.| 
2014ļand references therein). Near and mid-IR spectra were 


modelled by|Lisse et al.) (2012) to reveal large quantities of 
amorphous carbon and water ice-rich dust in the inner ring, 
as well as silicates. They interpret this as the recent break- 
up of a large centaur-like icy body due to a collision with an 
even larger silicate body, resulting in most of the current ob- 
served warm dust. Inwards perturbation of icy bodies from 
the outer cold belt by undetected planets as simulated by 
is supported by the large gap between 
the inner and outer disks, plus the subsequent detection of 
CO at ~ 20 au which could be due to inwardly evolving ex- 


ocomets undergoing CO sublimation (Marino et al.|2017). 


3.3. Late stages of stellar evolution 


3.3.1. Exocomets during the late stages of stellar evolu- 


tion 


Once a star evolves off the main-sequence, the large 
changes in luminosity in the Red Giant/Asympototic Gi- 
ant Branch (RG/AGB) stages, together with the mass-loss 
at the end of the AGB phase, will result in significant phys- 
ical and dynamical changes to any orbiting cometary bod- 
ies (Veras|2016). For example, during the Sun’s RG phase 
the water sublimation boundary will move from ~ 3 au 
to ~ 230 au, leading to cometary activity throughout the 


Trans-Neptunian region (Stern et al.|1990 2003). 
(2001) used the Submillimeter Wave As- 


tronomy Satellite (SWAS) to discover a large amount of 
water vapor around the carbon-rich AGB star IRC+10216. 
This was initially interpreted as originating by ongoing sub- 
limation of exocomets in that star’s Kuiper Belt. The in- 


ferred mass of cometary bodies was ~ 10g (Saavik Ford 
and Neufeld\2001). However, a later survey of eight AGB 


stars with Herschel showed a cometary 
source was unlikely, as the water line profiles matched that 
of the circumstellar CO, with widths as expected from the 
circumstellar expansion velocity. If the water originated 
from Kuiper belts, then at least some would be viewed near 
face-on, and the line widths of cometary-produced water 
would be much narrower. Instead, water formation appears 
to be either by shock-induced dissociation of CO, where 


the oxygen subsequently combines with hydrogen to pro- 
duce H20, or due to UV photodissociation of CO and SiO 
to again produce oxygen and water vapor 
(2016). Therefore, there currently exist no strong indicators 
of exocomets around stars at the RG/AGB stage. However, 
we know they are present due to observations of stars at the 
very end of their stellar evolution, White Dwarfs. 


3.3.2. Evidence for small bodies around white dwarfs 


There is now undisputed evidence that planetary sys- 
tems can (in some form) survive post-main sequence evo- 
lution. The origin of heavy element absorption lines in 
White Dwarf (WD) photospheres is now accepted as be- 
ing caused by pollution from infalling planetary mate- 
rial. The paradigm is that small bodies near the WD un- 
dergo tidal disruption to form a debris ring, where sub- 
sequent collisions create dust that is slowly removed via 
the Pointing-Robertson effect to “pollute” the stellar photo- 
sphere. Strong evidence for this is found in the association 
between WD infrared excesses due to circumstellar dust 
and atmospheric pollutants (see and references 
therein). 

The relatively short timescales for this material to sink 
below the observed photosphere means this process is cur- 
rently occurring where observed. The simplicity of WD 
photospheres results in high-precision elemental abun- 
dances being derived via model atmosphere fitting. There is 
no doubt that almost all small bodies providing this material 


were volatile-poor (Gdnsicke et al.2012} [Xu et al.|2019). 


However, two WDs have been identified as being polluted 
by water-rich bodies Raddi et al.{2015). 
This dominance of silicate-rich small bodies may not re- 
flect the original small-body population that survives into 
the WD phase. A theoretical study of water retention by 
found that larger water-rich 
bodies should survive the RGB and AGB stage of stellar 
evolution. They propose that the observed deficit is due to 
water being lost prior to tidal disruption or during the disk 
formation stage. If the WD is in a wide binary system, then 
dynamical evolution of exocomets from a Kuiper-belt via 
Kozai-Lidov resonances can result in WD pollution over a 
wide range of timescales of 108 — 101° yr 
(2017). Given that 25-50% of all WDs exhibit pollution by 
small bodies, exocomets may therefore be relatively com- 
mon around white dwarfs. 


4. DYNAMICAL EJECTION AND EVOLUTION OF 
COMETARY BODIES 


The discovery of 1I/‘Oumuamua and 2I/Borisov gives 
evidence to the hypothesis that over a star’s lifetime, part of 
its remnant small body population becomes unbound (e.g., 
Wyari2008). 
Thus ISO production seems to be a natural consequence of 
the existence of reservoirs of small bodies in planetary sys- 
tems. These small bodies can become unbound from their 
parent system in several ways (see Fig. [6}. In the following, 
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we will describe the various ISO production processes in 
chronological order of their occurrence during a star’s life- 
time. We give estimates on the order of magnitude level of 
the total mass of ISOs released by the individual processes. 
However, these values are highly sensitive to the assumed 
(debris) disk masses. 


4.1. 
4.1.1. ISO formation during disk dispersal 


In §2.2| we saw that disk masses decrease rapidly within 
the first 1-3 Myr. While initially, this mass loss will be in 
the form of small dust particles, later on, larger particles 
and eventually planetesimals will contribute to this mass 
loss. So far, the matter that becomes unbound during this 
early stage has received little attention. The main challenge 
is the significant uncertainties concerning the timescale on 
which planetesimals form. Therefore, it is unclear how 
many planetesimals are present in disks at this stage. Very 
young disks have masses comparable to the mass of their 
host star (mq ~ Ms). Considerable amounts of ISOs could, 
in principle, already be produced during this stage. The 
disk masses decrease to (mq ~ 0.01M/,) within just a few 
ten thousand years; if any planetesimals are already present 
during that phase, they could become released due to the 
lower potential caused by the gas loss. 

In recent years, several old disks ( > 10 Myr) have been 
discovered still containing considerable amounts of gas. 
These disks likely also contain a sizable planetesimal pop- 
ulation. When they eventually disperse their gas, the po- 
tential energy change due to the gas loss might be enough 
that a fraction of the planetesimals becomes unbound from 
the outer disk regions. Most proposed disk dispersal mech- 
anisms are relatively gentle. Thus the ISOs would leave 
their parent star at relatively low velocities (« 0.5 m/s). 
The only process that would potentially lead to asteroidal 
ISOs ejected at higher velocity are those involving disk jets. 
However, as jets only exist in the very early phases of disk 
development (<1 Myr) it is unlikely that planetesimals had 
already formed. 


Individual ejection mechanisms 


4.1.2. Cluster environment induced planetesimal release 


Planetesimals can also be released during another star’s 
close flybys (Hands et al.[2019). Here mostly planetesimals 
of the outer disk become unbound so that the ISOs produced 
by this process are usually icy objects. Close flybys happen 
most frequently during the first 10 Myr of a star’s life when 
the star is still part of a cluster (e.g. [Adams et al.\/2006). 
However, debris disks can also be affected, as their longer 
lifetimes can sometimes outbalance the lower frequency of 
such events during later stages. 

The frequency of close flybys varies enormously be- 
tween different types of stellar groups. In long-lived com- 
pact open clusters, stellar flybys are much more common 
than in the short-lived, more extended clusters typical for 
the current solar neighborhood. Based on typical disk trun- 
cation radii in diverse cluster environments, |Pfalzner et al.| 
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Kornmesser). 


find that clusters like the Orion nebula cluster are 
likely to generate the equivalent of 0.85 Mg of ISOs per 
star. In contrast, compact clusters like NGC 3603 can pro- 
duce up to 50 Mẹ of ISOs per star. Our solar system prob- 
ably created the equivalent of 2-3 Mẹ of ISOs by this pro- 
cess. In clusters, ISOs leave their parent system typically 
with velocities in the range of 0.5—2 km/s. 

External photo-evaporation and viscous spreading can, 
under certain circumstances, strongly influence the gas 
component of discs 
(2019). However, the decoupling of the gas and planetesi- 
mals means that these two processes release only a small 
number of ISOs. The velocity of these few ISOs should be 
much smaller than those released due to close flybys. 


4.1.3. Early scattering by planets 


During the phase when a planet is still accreting, plan- 
etesimals can receive gravitational kicks from planets. 
Such a perturbation can lead to two outcomes: scatter- 
ing and collisions. The ejection to collision rates (assum- 
ing constant density for the planet) are given as Re/Re = 


(M3 a2)(M2), where Mp and M. are masses of the star 
and the planet, and a is the planet’s semi-major axis 
(2018). At large distances, collisions become less likely, 
while scattering and, therefore, ejection becomes more ef- 
ficient. The reason is that the Hill radius, which controls 
scattering, expands linearly with the orbital radius, whereas 
the physical size is fixed. 

A planet can eject a planetesimal if the Safranov num- 
ber vesc /Vorb >> 1, where vesc is the escape velocity and 
Vorb the orbital velocity of the planetesimal 
[2018b). While the other giant planets all have similar 
Safronov numbers, it is actually the structure of the system 
itself that makes Jupiter responsible for the vast majority of 
planetesimals that have been ejected from the solar system 
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(Dones et al.}2015). 


Most investigations of this process assumed solar system- 


like environments (Raymond et al.||2020). The early so- 


lar system likely had a planetesimal disk with a mass of 


~ 5-65 Mg [Li er al |2022}. The 
estimates of how many planetesimals have been ejected 
range from 0.1 Mẹ (Raymond et al.|[2018b) to 20 Mg 
and even 40 Ma (Do et al.|/2018) 


per star. The differences in estimated total masses are pri- 
marily due to the unknown nature of the size distribution of 
ejected planetesimals. The ejected planetesimals from these 
regions would leave the systems with a typical velocity of 
4-8 km/s (Adams and Spergell2005), significantly higher 
than for the above mechanisms. The amount of planetesi- 
mals ejected by scattering depends on the structure of the 
planetary system. Systems of densely packed giant plan- 
ets that orbit close to a common plane are most efficient in 
ejecting planetesimals. 

In our solar system, the ejected planetesimals would 
have been primarily icy because ejection is much more ef- 
ficient outside the gas giants than closer in. Planetesimal 
ejection from the asteroid belt will be much less common. 
It is mainly caused by gravitational perturbations by Jupiter 
and Saturn and the mutual perturbations amongst the largest 
asteroids. It is currently an open question whether the early 


asteroid belt was high-mass and depleted (e.g. 
2001) or whether it was born low-mass and later on popu- 


lated by dynamical processes (Raymond and Izidoro|2017). 
If it was initially of high mass, ejections have depleted the 
asteroid belt by a factor of ~ 100 in number, the original 
asteroid belt being ~ 10-20 times more massive than today 
(O’Brien and Greenberg|2003). In this case, the total mass 
of ejected asteroids would amount to ~ 0.004 — 0.008 Mg. 

However, the system’s giant planets eject ISOs not only 


while still accreting gas (Raymond and Izidoro|2017} 


gies Zwart et al.}2018). In compact planet configurations, 
interactions between the planets themselves or the planets 


and the planetesimal disk can trigger instabilities. The plan- 
ets move positions throughout such a phase of planet-planet 
scattering. During such an event, the bulk of planetesimal 
disk would be ejected (Raymond et al.|2010). The structure 
of today’s Kuiper belt and the Oort cloud can be explained 
by this mechanism. 

Many exoplanetary systems are in compact configura- 
tions that might become unstable. Moreover, the solar sys- 
tem’s giant planets were likely originally on more compact, 
[2008). Originally it was proposed that such an instability 
happened in the solar system well after planet formation 
finished (500-700 Myr); however, a consensus is emerging 
that the instability happened early, no later than 100 Myr 
after the formation of the solar system (Zel/ner|2017} 
bidelli et al.|2018). 


4.1.4. Drifting from Oort clouds 


Even though the solar system’s Oort cloud formation still 
has many unknowns, we can also expect exo-Oort clouds 
to exist around other stars. Oort cloud bodies are icy exo- 


comets and possibly asteroidal objects 1999 
2016) and are only weakly bound to their par- 


ent system. Being situated at a significant distance from 
their parent star, they are subject to external influences, 
like Galactic tides and stellar flybys. Consequently, some 
planetesimals will always drift gently away from the star’s 
Oort cloud over the entire main-sequence lifetime of the star 


(Moro-Martin|2019b}|Portegies Zwart|2021). Thus their es- 
cape velocity is, on average, very low, probably < 0.1 m/s. 


(2018) estimates the Sun will lose 25-65% of 
its Oort cloud’s mass mainly due to stellar encounters. Such 
former Oort cloud members will contribute but amount to 


< 10% of the total ISOs population (Portegies Zwart|2021). 
4.1.5. White dwarf phase 


Finally, there is an increased planetesimal ejection rate 


towards the end of a star’s giant branch phase (Veras}2016 
[Veras et al.|[2020). During that phase, expansion of the 


stellar envelope can directly engulf closely orbiting plan- 
ets and tidally draw into the envelope planets which re- 
side beyond the maximum extent of the stellar envelope. 
While the planetesimal disks are mainly unaffected by en- 
gulfment, stellar mass loss plays a role. During the post- 
main sequence evolution, all objects, planets and planetes- 
imals alike, move outward due to orbital expansion caused 


by mass loss (Villaver et al.||2014). In the solar system, 


planetesimals within ~10° au would all double their semi- 
major axes and therefore be less strongly 
bound to their host star. The ISOs produced during the gi- 
ant branch phase gently drift away from their parent star 
at very low relative velocities of 0.1 — 0.2 km/s 
[et al./2021a). In addition, mass loss also destabilizes the 


planetary system, resulting in the gravitational scattering of 


14 


planetoids by massive perturbers. (2018) suggests 


tidal disruption events of (initially bound) planetary objects 
by white dwarfs as an additional ISO production mecha- 
nism. Some of these objects are scattered towards the white 
dwarf on almost radial orbits and get tidally shredded apart. 
Some of the fragments that are produced are ejected into 
interstellar space. 

The highest mass white dwarf progenitors would yield 
the greatest giant branch excitation and and are most effi- 
cient in the ejection of planetesimals (Veras et al.\/2020). 
This study points out that most stars in the Milky Way are 
less massive than the Sun, and concludes that the produc- 
tion of ISOs from within 40 — 1000 au of the evolved star is 
insignificant when compared to the number created during 
early stellar and planetary evolution. 


4.1.6. Binary scenarios 


In we noted that binary and multiple stars are very 
common in the Galaxy. Moreover, binaries are also more 
efficient scatterers, providing more substantial dynamical 
perturbations than even a gas giant planet (Raymond and 
cesses discussed before for single stars could also take place 
in binary systems. So far, more or less exclusively scatter- 
ing processes have been considered in detail. Close binary 
stars have well-defined dynamical stability limits, and any 
objects entering within a critical orbital radius are desta- 
bilized. assume that planetesimals 
form beyond this stability limit and drift inwards across 
the limit due to aerodynamic gas drag. Their simulations 
show that all planetesimals that are drifting inside the sta- 
bility limit are ejected. Especially for close binaries, the 
snow line is often exterior to the stability limit such that 
a considerable portion of their ejected planetesimals may 
be refractory. However, this mechanism is only efficient if 
a significant fraction of the solid disk mass drifts interior 
to the binary stability limit. For the aerodynamic forces 
to work efficiently the planetesimals have to be relatively 
small (r < 1 km). Only then can they enter the dynamically 
unstable zone close to the binary, rather than piling up at 
the pressure maximum in the gaseous circumbinary disk. It 
has been proposed that 11/‘Oumuamua could be a fragment 
of a planet that formed in a binary system and was tidall 
disrupted after passing too close to one of the stars 
2018). However, tidal disruption is a rare event. Neverthe- 
less, tight binary systems can eject an amount of rocky ma- 
terial comparable to the predominantly icy material thrown 
out by single- and wide binary-star systems. It is estimated 
that the mass of ISOs ejected from binary systems is at least 


the same as by planet scattering (Jackson et al.|2018) and 


the ISOs leave their parent binary system with ~ 6.2 km/s 


(Adams and Spergel\2005). 


4.1.7. Relative importance of the different ISO formation 
mechanisms 


All mechanisms mentioned in this section produce ISOs. 
The question is which of them is/are the dominant ISO pro- 
duction mechanism(s). All the processes discussed above 
produce of the order of 0.1 — 30 Mga per star. However, 
given the strong dependence on the assumed mass of the 
planetesimal reservoir, these numbers fail to give a clear 
winner. Thus different criteria have to be employed. The 
various models make different predictions concerning the 
properties of the ISOs. As soon as a large enough sample of 
ISOs is known, these model predictions can be used to iden- 
tify the dominant ISO production mechanism. Differences 
are mainly in terms of the composition and the velocity of 
the ISOs. 

Most mechanisms lead predominantly to icy ISOs that 
formed outside the snowline in their parent systems. How- 
ever, 1J/‘Oumuamua did not show the cometary coma mor- 
phology one would expect from an icy ISO. Therefore, 
many mechanisms have been tested concerning their ability 
to produce refractory ISOs. In general, scattering processes 
tend to produce more asteroidal ISOs than stellar flybys, 
drifting from the Oort cloud and the processes at the end 
of the main sequence. However, even scattering processes 
produce mainly icy ISOs. 

There are significant differences in the ejection velocity 
of the various processes (Pfalzner et al.|202 1a). The slowest 
ISOs are those that drift away from the Oort cloud or are 
shed during the stellar post-main-sequence phase (~ 0.1- 
0.2 km/s). Stellar flybys lead to ejection velocities in the 
range of 0.5-2 km/s; even higher velocities are achieved by 
planet scattering (~ 4-8 km/s), only exceeded by processes 
where two or more stars kick out the planetesimals together. 

Most processes discussed have a connection between ve- 
locity and composition. Due to their closeness to the star, 
the refractory planetesimals require a higher velocity to be- 
come unbound than the icy planetesimals (unless, of course 
the rocky planetesimals resided in the Oort cloud of the host 
star). As soon as a statistically significant sample of ISOs is 
discovered, the combined information of their observed ve- 
locities and composition might help constrain the dominant 
production process. However, the observed ISO velocity is 
only indirectly linked to the ejection velocity, and other ef- 
fects such as dynamical heating for older ISOs may make 
such studies difficult. The complete ISO velocity distribu- 
tion contains multiple components reflecting various ejec- 
tion speeds and the parent system’s different ages. In fu- 
ture, disentangling the different components will be one of 
the significant challenges. 


4.2. Resulting ISO population 


Most of the above described planetesimal ejection mech- 
anisms produce considerably more icy ISOs than rocky 
ones. The ejection of volatile-free asteroids is possible, but 
they are thought to be a relatively small fraction of planetes- 
imals that were ejected. However, this impression might be 
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partly caused by our heliocentric view. In the solar system 
giants planets all orbit beyond the “snowline”. Therefore, 
the majority of bodies within <2.5 au on unstable orbits 


end up colliding with the Sun (Gladman et al.|1997}{Minton| 


and Malhotra|2010). In the solar system system there is in- 
deed a much greater supply of icy planetesimals, however, 
this may not apply elsewhere iCutlo018), If we look at 
the known exoplanet systems, many of them have relatively 
massive planets much further in than in our solar system. 
The ejection of rocky planetesimals would be more effi- 
cient where there are close-in massive planets around high- 
mass stars. However, even then it appears unlikely that the 
rocky planetesimals could dominate the galactic population 
of scattered ISOs. The reason is that the large reservoir of 
volatile planetesimals beyond the snowline is always much 
larger than that of the rocky planetesimals in the inner sys- 
tems. 

Naturally, the size distribution of ISOs is connected to 
the size distributions of the reservoir of planetesimals the 
ISOs originate from. Unfortunately, there is still a large 
uncertainty of the size distribution of planetesimals. Often 
a mass function of the form Nrso, x m`? is assumed. 
This means that the number of ejected particles of a certain 
mass N(m) is then given by 

_2-p_ Mrso 


) p- 1m-imip” 


(Adams and Spergel|2005), where Mrso is the total mass 


of ISOs and m,,, the largest mass possible. Sometimes p = 
5/3 is chosen (Adams and Spergel|2005). However, starting 
from a Dohnanyi distribution and from a distribution ob- 
tained by simulations of the streaming instability (p=1.6), 
argue that the underlying masses 
are off by 2-4 orders of magnitude. Observational con- 
straints of the small body population of the solar system 
can also provide some constraints on the initial size dis- 
tribution of ISOs. Beyond Neptune, the small body popu- 
lation is thought be the least collisionally evolved 
let al.[2022). For these planetesimals a turnover in the size- 


number distribution is observed around a radius of 50 — 70 
Bernstein et al. [2004 , which likely 
is of primordial origin (Vitense et al./2012). Therefore, one 
can expect that size distribution of the ISOs before expul- 
sion was similar. 

The size distribution of the planetesimals before expul- 
sion sets the stage for the ISO size distribution. How- 
ever, both distributions are not necessarily identical. First, 
there might be selection effects concerning the planetesi- 
mals mass (and therefore size) in some of the ISO pro- 
duction mechanisms. For example, in three body scatter- 
ing processes, ejection is more likely the lower the mass of 
the body. Second, the ejection processes themselves can al- 
ter ISOs, such that the properties of the ISOs differ from 
those of the planetesimals in the parent system. While ISOs 
that are relatively gently released from their parent system 
during the Oort cloud phase will remain unaltered, ISOs 
ejected in a more violent fashion might be disrupted during 
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Fig. 7.— Top: Stellar age distributions in the Milky Way, 
adapted from [Gallart et al.|(2019). Bottom: Age distribu- 
tion of planet hosting stars as listed in the public database 
at exoplanets.eu. The observational bias towards surveying 
Sun-like stars is clearly evident. 


ejection (Raymond et al.\|2020). For example, planetesi- 


mals gravitationally scattered in very close encounters with 
a giant planet can be subjected to tidal disruption (e. g.{As-| 
haug and Benz) 1996). In particular, scattering events that 


include a close passage to the parent star can lead to the de- 


cline or loss of cometary activity (e.g. 
1997), and sublimation-driven activity may also flatten fu- 


ture ISOs’ shapes 
(2021). Third, ISOs might be affected to different degrees by 
their journey through interstellar space (see $4.4). 

By now, it should have become clear that planet forma- 
tion and the production of ISOs are intrinsically linked pro- 
cesses. Ever since planets formed, a population of plan- 
etesimals were produced at the same time. A considerable 
portion of these planetesimals are released sooner or later 
from their parent star(s) and became ISOs. This means that 
the interstellar medium is steadily enriched with newly re- 
leased ISOs (Pfalzner and Bannister|2019). Just like the 
Galaxy’s stellar metallicity increases over time so does the 
density of interstellar objects. 

The ISO age distribution should be directly linked to the 
age distribution of the planets and stars. The stellar age 
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distributions within the Milky Way are shown in Fig. [7| It 
demonstrates that star formation peaked at the early ages 
of the Galaxy. However, it is unclear whether the first 
stars produced ISOs at the same rate as today. It is known 
that planet formation is much more efficient around high- 
metallicity stars. ISO formation is linked to planet forma- 
tion. Suppose we take the age distribution of the planet host 
stars (Fig. [7|bottom) as a proxy for the planetesimal produc- 
tion distribution, this differs considerably from the general 
stellar age distribution in the Milky Way. The age distribu- 
tion of planet host stars shows a clear maximum at approxi- 
mately 4 Gyr. The age distribution of ISOs might shift even 
more to younger ages, as some portion of the ISOs may dis- 
solve during their interstellar journey since the frequency of 
volatile loss and extinction is far higher for ejected planetes- 


imals than for surviving ones (Raymond et al.|2020). 


4.3. 


A planetesimal that has been ejected from its host star 
and has become an ISO does not necessarily remain an ISO 
forever. ISOs may be recaptured. Recapture is most likely if 
the planetesimal has been ejected from a host star which is 
a member of a star cluster. This situation is most likely for 
young host stars because most stars are born in a star cluster 
environment (Lada and Lada\2003). However, some stars 
remain for hundreds of Myr in open clusters and for several 
Gyr in globular clusters. In all these situations, recapture is 
a real option. 

In such an environment, the ISO might be very quickly 
(< 1 Myr) recaptured by a different member of the same 
star cluster and lose its status of an ISO (Hands et al.[2019). 
The likelihood of the ISO being recaptured rises with the 
cluster mass as the escape speed from the cluster increases 
with cluster mass. Equally, ISOs escaping from stars in the 
cluster centre are more likely to be recaptured. For ISOs 
that are not in a star cluster environment, the likelihood 
of recapture is much lower. In any case, the velocity of 
the ejected ISO is the key parameter for being recaptured. 
The lower the ISO’s velocity the higher the probability that 
it is recaptured. This correlation is valid in every capture 
environment, whether it is another star, a disc surround- 


ing a young star or a molecular cloud (Grishin et al.)2019 
Pfalzner et al./2021a 2022). 


4.4. Physical processing in the ISM 


ISO timescales 


ISOs are often viewed as travelling through interstel- 
lar space in calm cryogenic conditions. Therefore, ISOs 
are expected to be pristine samples of the planetesimals of 
other planets. As planetesimals are leftovers from the planet 
formation process, ISOs should therefore give us direct in- 
formation about planet formation elsewhere. Generally, an 
ISO’s journey through interstellar space might be not as un- 
eventful as often imagined. The environmental effect on 
ISOs can be expected to be similar to that of Oort cloud 
comets. Just like the comets residing in the Oort cloud 


(Stern|2003), a variety of thermal, collisional, radiation, and 
ISM processes might affect ISOs during their long journey 
(Stern|1990). The ISOs’ surface is subject to these environ- 
mental influences, which can potentially lead to the erosion 
or even complete destruction of an ISO. How destructive 
these processes are depends on the composition of the ISO, 
with a rocky ISO being much less affected than, for exam- 
ple, a hydrogen-rich ISO. 

During their journey, ISOs experience cosmic mi- 
crowave radiation and radiation from the stars. The tem- 
perature increase by the cosmic microwave radiation on 
the surface of the ISO is probably negligible. In contrast, 
the effect of stellar radiation depends strongly on the ISO’s 
individual journey. Any ISOs passing close to a star can 
be expected to be modified considerably by such an en- 
counter. The prime effect is the reduction or complete loss 
of volatiles. Such an event is expected to be rare. However, 
parsec-range and closer encounters with highly luminous O 
and supergiant stars can heat ISO surfaces to temperatures 
capable of removing the most volatile ices, such as neon or 
oxygen (Stern|2003). 

ISOs are also subjected to the effects of cosmic rays and 
stellar energetic particles. Both have a broad spectrum of 
energies and interact with the ISO surface and subsurface. 
Cosmic rays are the primary source of space weathering for 
the comets in the Oort Cloud and, therefore, likely also for 
ISOs. While low-energy particles interact only with the 
ISO surface, the most energetic ones deposit a significant 
amount of energy down to tens of meters. These processes 
can modify the isotopic ratios in cometary ices and cre- 
ate secondary compounds through radiolysis 
(2020). The penetration depth of cosmic rays is relatively 
small in comparison to the ISOs of size > 100 m that are 
likely to be detected by telescopes, see §5| Even for the 
extreme case of pure hydrogen ISOs, it would be < 10 m 
(Hoang and Loeb\|2020), for rocky material considerably 
shallower. Thus cosmic rays immediate effect on ISOs is 
small, but over time it can erode an ISO gradually from the 
surface or affect the composition of the upper layers (see 
97.3.2). 

The interstellar medium contains large amounts of gas 
and dust. The effect of this gas and dust on ISOs depends on 
the gas and dust densities and the relative velocity of the in- 
dividual particles. Gas and dust densities are relatively low 
in the interstellar medium itself (ng ~ 10 cm~®) and are 
considerably higher in molecular clouds (ng ~ 104 cm~?) 
and even higher when passing through protoplanetary disks 
of newly forming stars (ng œ~ 10577 cm~?). 

Collisions of ISOs with the ambient gas at high speeds 
can heat the frontal area, possibly resulting in transient 
evaporation. Such a situation occurs when the ISO passes 
through a molecular cloud. Each particle directly impacting 
on the surface of an ISO delivers an energy Ep = mpv? /2. 
Therefore, ISOs are most impacted by high-velocity parti- 
cles. At least for Oort cloud comets, simulations show that 
dust particles could significantly erode the cometary sur- 


face in the range of several meters 1990 
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1993), preferentially removing the sub-micron grains. 


However, erosion by ISM grains is a complex process de- 
pending on the ISM grains composition and structure. 

While moving through a relatively massive molecular 
cloud forming a large star cluster, there is a non-negligible 
probability that the most massive stars explode as a super- 
nova. Therefore, an ISO might be subjected to the explo- 
sion products from such a supernova. Despite supernovae 
explosions being brief (~ 0.1 yr), they are extremely lumi- 
nous L = 1092.4, and therefore can heat ISOs from con- 
siderable distances. Their intense but shorter thermal pulses 
could propagate 0.1-2 m into ISO surfaces (Stern|2003). 

These processes complicate the interpretation of ISO ob- 
servations; the surface layer of ISOs might be modified over 
time. In summary, all these processes erode the surface of 
ISOs to different degrees. However, as the distinctive prop- 
erties of the ISOs in terms of their size and composition are 
not well constrained, it remains unclear whether their sur- 
face is affected or whether they decrease in size with age 
and how many small ISOs become destroyed. 


5. ISO SEARCHES AND DISCOVERY 


5.1. 


Upon discovery, ISOs will appear as either comets or as- 
teroids depending on their level of activity. Hence the single 
most important discriminant is showing the orbit is hyper- 
bolic, with e > 1 anda < 0 at the 3-ø level or higher. 
Dynamically, this is equivalent of saying that their veloc- 
ity at an infinite distance vso > 0 (also known as the hy- 
perbolic excess velocity). This requires accurate astromet- 
ric measurements of the object over a period of time, and 
hence sensitive wide-field surveys that allow early detec- 
tion. For ISOs close by the Earth it may take only a few 
days of observations before the orbital arc is large enough 
to show this unambiguously; more distant ISOs may take 
weeks (see $6). 

It should be noted that many comets are discovered each 
year on apparently hyperbolic orbits. These are in real- 
ity dynamically new comets from the Oort cloud, but are 
listed as hyperbolic due to a number of reasons. First, an 
orbit that is near-parabolic in the solar system barycentric 
reference frame may appear to be hyperbolic in the he- 
liocentric reference frame, in which most orbital elements 
are published. For example, the dynamically new comet 
C/2021 A9 (PANSTARRS) has an osculating heliocentric 
orbital eccentricity of e = 1.004. Transforming to the cor- 
rect barycentric frame gives e = 1.0004, much closer to 
a parabolic orbit. Secondly, an originally weakly bound 
object may become unbound due to gravitational pertur- 
bations or non-gravitational forces due to outgassing, see 
[Królikowska and Dybczyrski (2010). 

To illustrate this, using data from the JPL Small Bodies 
Database ( https://ssd.jpl.nasa.gov/tools/sbdb_query.html ), 
17 Long-Period Comets (LPCs) were discovered in 2020- 
2021 with osculating orbital eccentricity e > 1.0. Yet all 
eccentricities were so close to unity that either these comets 
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have e < 1.0 within the orbital uncertainties and/or within a 
barycentric reference frame, or they would have possessed 
an originally parabolic/elliptical orbit before entering the 
our planetary system. Understanding these factors is crucial 
for confirming the status of suspected ISOs with e ~ 1. 


5.2. 


There were many papers that made predictions about 
the space density of ISOs based on the non-detection in 


surveys. Some very early work was reported by 
(1972), (1975); [Sekaninal (1976): 
(1987); Valtonen and Zheng (1990); [Zheng et al.| (1990). 
The advent of wide-field CCD-based sky surveys in the 
1990s allowed quantitative upper limits to be assessed. 


(2005) derived a 95% upper confidence limit of 
4.5 x 1074 au~? using the LINEAR survey. Just before 


the discovery of 1J/‘Oumuamua, 
numerically integrated ISOs and linked them to the non- 
detection by the combined major sky surveys operational at 
that time; the Catalina Sky Survey, the Mt. Lemmon Sur- 
vey and the Pan-STARRS project. Although they assumed 
isotropic approach trajectories which disregarded the local 
standard of rest distribution of stellar velocities, they pre- 
sciently derived upper limits for both inactive and active 
(cometary) ISOs. They reported 90% confidence upper lim- 
its for 1 km diameter ISOs of < 2.4 x 107? au~? and 
< 1.4 x 1074 au~? respectively. 

The important thing to get across here is why was the 
first ISO discovered now? To discover these fast moving 
objects you need to survey the whole sky quickly it to very 
faint limiting magnitudes, and we have only had surveys 
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capable of doing this in the last couple of decades. [Heinze] 
shows that current near Earth object (NEO) 
surveys are poorer than expected at finding fast moving ob- 
jects. However, it is expected that future surveys will do 
much better (see 97.3.1). 


6. OBSERVED CHARACTERISTICS OF ISOs 
6.1. 


Some of the most important questions about the first ISO 
include: (1) Where did it come from? and (2) What is it 
made of? As discussed below, we still don’t know the an- 
swers to these questions. 1I/‘Oumuamua was easily observ- 
able from the ground for a little over a week, but as it moved 
away from the Earth it faded quickly, and large telescopes 
were able to observe for only about 1 month. The last Hub- 
ble Space Telescope observations were made in Jan. 2018. 
It is remarkable that we know as much about this object as 
we do, because all of the large telescope time had to be se- 
cured through Director’s requests. In total approximately 
100 hrs on 2.5-10-m ground-based telescopes were devoted 
to characterizing this exceptional object. To date, over 200 
refereed papers have been written on both interstellar ob- 
jects, and nearly 450 papers including non-referred mate- 
rial. It is not practical to cite everything, so key papers have 
been highlighted, and we have included some reviews of the 
field. 


1I/‘Oumuamua 


6.1.1. Discovery of 1I/‘Oumuamua 


On 2017 Oct. 19 the Pan-STARRS survey found an ob- 
ject, designated P10Ee5V, moving quickly with respect to 


Fig. 8.— Images of the two ISOs passing through the inner solar system show the dramatic difference in their characteris- 
tics. Top row: 1I/‘Oumuamua: (a) Pan-STARRS discovery image (top) taken on 2017 Oct. 19, with follow up images from 
the CFHT near min/max brightness on Oct. 27; (b) Gemini 8-m color composite image made from 192 images (1.6 hrs) 
showing no hint of a dust coma (Gemini Observatory/AURA/NSF); (c) artist’s depiction of two possible nucleus shapes 
based on the large light curve range (Credit: ESO, M. Kornmesser; William Hartmann) and artist’s view of the ISO after 
the discovery of non-gravitational acceleration (ESA/Hubble, NASA, ESO, M. Kornmesser). Bottom row: 2I/Borisov: (d) 
CFHT image taken 10 days post-discovery (2019 Sep. 10); (e) HST image on 2019 Oct. 12 (NASA/ESA)); (£) radially 
normalized ratio HST images median averaged over an orbit from 2019 Oct. 12, Dec 24. (00:31 UT), and Dec 24. (02:21 
UT) showing outgassing jets, and 2020 Jul. 6 showing the split nucleus (NASA/ESA), processed by H. Boehnhardt. 


the stars (see Fig. Bb during the normal near earth ob- 
ject survey. Rob Weryk then found pre-discovery images 
in Pan-STARRS data from Oct. 18. Follow up data ob- 
tained by the 1m ESA ground station on Oct. 20 was re- 
jected by the Minor Planet Center because the data im- 
plied a large eccentricity, and P10Ee5V was classified as 
an Earth-orbit crossing asteroid. Data from Oct. 20 ob- 
tained by the Catalina Sky Survey suggested that the object 
should be classified as a short period comet. However, ob- 
servations obtained from CFHT on Oct. 22 showed that 
the orbit was hyperbolic with an eccentricity of 1.188. On 
Oct. 24 the object was designated C/2017 U1 (MPEC 2017- 
U181). This was corrected on Oct. 26 after deep images 
from the CFHT from Oct. 22 showed no coma. The object 
was named A/2017 U2 (MPEC 2017-U183). 

Within a week of discovery a request was made to the 
Hawai’ian cultural group that advised Maunakea observa- 
tories management group to propose a name for the new 
object. They proposed the name ‘Oumuamua, meaning “a 
messenger from afar arriving first” or a scout or messenger 
sent from our distant past to reach out to us or build con- 
nections with us. In an extraordinarily fast effort, the IAU 
approved this on 2017 Nov. 6, and the new name became 
1I/“Oumuamua. This became the foundation of a project 
called A Hua He Inoa which blends traditional indigenous 
practices into the official naming of astronomical discover- 


ies (Kimura et al.[2019}|Witze|2019). 


6.1.2. Nuclear Characteristics 


One of the most straight forward measurements to make 
is the brightness, and from this get an estimate of the ob- 
ject’s radius for an assumed albedo: 

pray = 2.235 x 1072r2 A719 4(mo-™)1 199-4062) (2) 
where r, and A are the helio- and geocentric distances 
[au] and mọ and m are the apparent magnitudes of the 
sun and comet (Russell|[19 16). It was apparent from the 
earliest imaging observations of 1I1/ʻOumuamua that it had 
a very large rotational light curve range so various esti- 
mates of the size were reported, even for the same assumed 
albedo. Combining several nights worth of observations 
from several observatories yields Hy = 22.4 which gives 
an average radius of 0.11 km assuming a typical cometary 
albedo of 0.04 (Meech et al.|2017b). Spitzer observations of 
11/“Oumuamua, in principle could have provided measure- 
ment of both the nucleus size and albedo, however because 
of strict solar avoidance angles these observations could not 
be made until late 2017 Nov. and only upper limits on the 
flux were obtained (Trilling et al./2018). Thermal models 
were used to estimate the corresponding radius and albedo, 
with a preference for an albedo of 0.1 and a radius consis- 
tent with the previous estimate. 

The rotational light curve shown in Fig. has a bright- 
ness range of Am ~ 3 mag, implying an axis ratio for an 
oblate sphere of ¢ = 10°4(4™ = 15, However, this does 
not take into account the phase angle (which can make the 
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object look more elongated than it is (see Fig. BE). It also 
does not take into account that the rotation pole position is 
completely unknown; if the pole was more closely aligned 
to the ecliptic, then the object could appear even more elon- 
gated that the light curve suggests. The consensus is that 


the ratio is likely a/b = 6 (‘Oumuamua ISSI Team et al. 
2019; |Mashchenko\|2019). Solar system objects typically 


do not have axis ratios this large, and the cause of the un- 
usual shape of 11/‘Oumuamua remains one of the enduring 
mysteries (see §7.3.2). 

Attempts to find a rotation period from the data pro- 
duced a variety of periods, near 8 hrs, all differing depend- 
ing on the length of the data set. The most comprehensive 
model concluded that 1J/‘Oumuamua is in a complex ro- 
tation state, rotating around its shortest axis with a period 
of 8.67+0.34 hours with a period of rotation around the 
long axis of 54.48 hours (Belton et al.2018). The damp- 
ing timescale for a body this small is long enough that an 
excited rotation state can be preserved from the time of ejec- 
tion from the host star. The shape of 11/‘Oumuamua as in- 
terpreted from the rotational light curve can change signifi- 
cantly, depending on the rotation state (see Fig. [8f). 

used the shape and rotation of 
1J/‘Oumuamua to place some constraints on the strength 
and density of the ISO, finding it was more likely to have a 
density typical of asteroids. 


6.1.3. Constraints on composition 


Several groups obtained spectra for 11/“Oumuamua, both 
from the ground and from Spitzer, but there was no evidence 
of outgassing (see 
for a summary). Many groups used either spectra or fil- 
ter photometry to estimate the spectral reflectivity of the 
surface. The spectral reflectivity, S, was found to increase 
with wavelength and ranged between S = (7-23)+3%/1000 
A (ewiri and Seligman|2022}. This measurement was chal- 
lenging because the rotational signature had to be removed. 
Regardless of the value, the surface of 1I/‘Oumuamua was 
red, typical of organic-rich comet surfaces, but also consis- 
tent with iron-rich minerals and space weathered surfaces. 

Deep stacks of images from many nights of data from 
ground based and space-based optical images showed no 
dust at all to a limit of there being < 10~° kg/s of micron- 


sized grains being produced (Meech et al.|/2017b 
[et al.|[2017). In contrast, the dust production upper limit 


from Spitzer observations for 10 um grains was < 9 kg/s 


(Trilling et al.\2018). 


6.1.4. Trajectory and potential origins 


One of the primary goals for the HST time awarded for 
11/‘Oumuamua was to extend the astrometric arc length to 
be able to determine its orbit with sufficient precision to 
trace its path backwards and determine its home star sys- 
tem. The chances of being able to trace the trajectory back- 
wards to find its parent star were low because gravitational 
scattering from stellar encounters would limit the search to 


2017 Oct 27 
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Fig. 9.— Rotational light curve of 1I/‘Oumuamua compiled from data taken from various observatories in 2017 October. 


the past few tens of millions of years (Zhang|2018). 


As 1I/‘Oumuamua faded, astrometry was obtained with 
the VLT in Chile and with HST, with the last observations 
from HST being obtained on 2018 Jan. 2 when the ISO was 
fainter than V ~ 27 at 2.9 au. An analysis of the com- 
bined ground and HST observations showed that the or- 
bit could not be well fit with a gravity only solution. By 
adding a radial non-gravitational acceleration term directed 
radially away from the Sun with a a dependence, the or- 
bit was well fit (Micheli et al.|/2018). The authors ruled 
out a number of other more hypotheses, e.g. Yarkovsky 
effect, frictional drag forces, impulsive velocity changes, 
photocenter offset from a binary or non-uniform albedo, in- 
teraction with the solar wind or radiation pressure. Most 
were rejected because the acceleration produced by these 
mechanisms was either too small, in the wrong direction, 
wouldn’t match the continuous nature, or implied a non- 
physical bulk density for 11/‘Oumuamua. By assuming that 
the same non-gravitational acceleration was operating on 
the pre-perihelion trajectory, an attempt was made to find 


the originating star system for 1I/“Oumuamua. |Dybczynski 
and Królikowska| (2018) attempted a search for the parent 


star taking advantage of the Gaia satellite data release 1 
(DR1) of stellar positions, but did not find an obvious par- 
ent star for 11/‘Oumuamua. 

The Gaia DR2 data release provided the necessary as- 
trometric position and velocity information that was needed 
to be able to try to trace the path of 11/‘Oumuamua back to 
its home star using its Keplerian orbit 
2018). Plausible home stars would be those which had 
been within 0.5 pc of 1]/‘Oumuamua’s trajectory (the ap- 
proximate size of the Sun’s Oort cloud), and which had low 
encounter (ejection) velocities < 10 km/s (see §4.1). An 
initial selection of stars which passed within 10 pc of the 
trajectory was selected and the candidates and the path of 
11/“Oumuamua were integrated through a smooth Galactic 
potential. This resulted in four potential “home” systems 
with encounter distances between 0.6-1.6 pc and encounter 
velocities between 10.7-24.7 km/s, all with ejection times 
less than 6.3 Myr ago. Unfortunately, all of these systems 
had ejection speeds much higher than expected for ISOs. 
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6.1.5. Cause of the non-gravitational acceleration 


In order to explain the non-gravitational acceleration of 
11/‘Oumuamua in the absence of any apparent activity (dust 
coma or gas), |Micheli et al.|(2018) estimated the mass loss 
needed to accelerate the nucleus to the observed value for a 
range of comet and asteroid densities. They found mass loss 
requirements between 0.7-140 kg/s, but adopted 10 kg/s as 
the best estimate. A thermal model matching the accelera- 
tion would require sublimation of water plus a more volatile 
species (such as CO% or CO). With a high CO/H20 ratio 
the model could produce enough mass loss to within a fac- 
tor of 2-3 of what was needed to accelerate 1I/‘Oumuamua. 
found 3¢ upper limits for the dust 
production of 9 kg/s (for 10 um grains), Q(CO2) = 9 x 
107? mol s~! and CO 1.1 x 1024 mol s~! (after correc- 
tion). None of the optical spectra were sensitive enough 
to detect gas with such a low CO or H20 production rate. 
Having consistency within a factor of a few between the 
model and observations suggests that volatile outgassing is 
a plausible scenario for the acceleration of 1I/‘Oumuamua. 
postulated that any comet passing through a 
molecular cloud would have the small grains eroded from 
the surface, and this could explain the lack of dust in the 
presence of outgassing. 

Other more exotic volatile sublimation scenarios to ex- 
plain the non-gravitational acceleration are discussed in 
7.3.2] Finally, investigated the 
spin dynamics of 1]/‘Oumuamua with outgassing consis- 
tent with the non-gravitational acceleration and found that 
this need not cause a spin up of the nucleus. They were able 
to reproduce the observed light curve with a model of CO 
outgassing but not with water-ice jets. 


6.2. 21/Borisov 
6.2.1. 


The second ISO was discovered on 2019 Aug. 30 by 
Gennadiy Borisov, an engineer at the Crimean Astronom- 
ical Station and amateur astronomer, using a 0.65m f/1.5 
astrograph (Borisov and Shustov|2021). Borisov reported it 
as a potential new comet with a compact 7 arcsec diame- 
ter coma, and detected a short 15 arcsec tail two days later. 
It was placed on the Minor Planet Center Potential Comet 
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Confirmation Page as object gb00234, and by 2019 Sep. 11 
sufficient additional astrometry had been reported to give a 
cometary designation of C/2019 Q4 (Borisov) and an initial 
orbit with an eccentricity e = 3.08 (MPEC 2019-R106). 
On 2019 Sep. 24 the permanent designation of 2I/Borisov 
was announced (MPEC 2019-S72). As further observations 
were reported, by mid-October the derived orbit had stabi- 
lized with e = 3.354, q = 2.006 au and ¿i = 44.05°, giving 
discovery circumstances of r, = 2.99 au and A = 3.72 au. 
This allowed [Fe et al.| (2020) to identify precovery images 
in archival survey data dating back to December 2018 when 
it was at a heliocentric distance of 7.8 au. 


6.2.2. Trajectory and potential origins 


The significantly longer observation arc for 2I/Borisov 
allowed the measurement of standard cometary non- 
gravitational forces, implying that searches for progenitor 
systems might suffer less uncertainty than for 1]/“Oumuamua. 


However, |Hallatt and Wiegert (2020) did not identify any 


systems with the golden combination of small encounter 
distance and low encounter velocity. 
investigated a range of non-gravitational force mod- 
els and found that 21/Borisov was ~ 0.053 — 0.091 pc from 
the MOV star Ross 573 910 kyr ago. They also demon- 
strated that such a close encounter is unlikely over this short 
time period. However the predicted encounter velocity was 
22.6 km/s, significantly higher than expected for ejection of 


ISOs, even assuming giant planet interactions (Bailer-Jones| 
et al.|2018). 


Hence like 1]/“Oumuamua, the home system of 21/Borisov 
remains unidentified. It is likely that unless the ejection was 
from a nearby star very recently, that tracing the path to the 
home system is not possible. As soon as an ISO passes 
through a molecular cloud, back-tracking no longer possi- 
ble. 


6.2.3. Nuclear Characteristics 


Due to its activity, there were no direct detections of the 
nucleus of 2I/Borisov. In their archival analysis, [Ye et al] 
did not detect 2I/Borisov in deep co-added images 
from November 2018. This gave a strong upper limit to the 
nucleus radius of ry < 7 km assuming a geometric albedo 
of 0.04. Rigorous imaging constraints came from HST ob- 
servations by Jewitt et al.|(2020a), who measured a strong 
upper limit of Hy > 16.60 + 0.04 assuming a nominal 
cometary phase function of 0.04 magnitudes/degree. Fit- 
ting the coma surface brightness distribution gave showed 
no significant excess from a central nuclear psf, giving an 
improved upper limit to the nucleus radius of r < 0.5 
km, or Hy > 19.0. calculated the non- 
gravitational parameters for a range of heliocentric distance 
laws, and combined them with outgassing measurements 
(see below) and an assumed nuclear density of p = 0.5 
gm cm™? to derive r < 0.4 km. used 
measured water production rates to derive a minimum ra- 
dius of r > 0.4 km assuming the entire sunward surface of 
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the nucleus was active. 

Several authors used apparent jet structures visible from 
HST imaging of the inner coma to attempt to constrain 
the spin axis. (2020) found it difficult to ex- 
plain the overall dust coma anisotropy in terms of a sin- 
gle isolated emission source, and preferred a model of gen- 
eral mass loss from the afternoon nuclear surface, giving a 
pole orientation of (a, ô) = (205°, +52°). From the same 
data, interpreted the appearance of 
two jet/fan structures as implying localized emission near 
the rotation equator, and derived a spin pole direction of 
(260°, —35°). aA eee a data with 
additional HST imaging, (2020) interpreted 
the positional evolution of the jets as implying a near-polar 
source with a direction (322°, +37°). Although all three 
studies reported uncertainties of at least +10°, the disparity 
between these results led to con- 
clude that it was not possible to sufficiently constrain the 
pole orientation with the extant data. The latter study also 
reported a tentative detection of a 5.3 hr periodicity in the 
light curve, but the low amplitude of ~ 0.05 magnitudes im- 
plies a low significance. The lack of precise determinations 
of these nuclear properties are all unremarkable when com- 
pared to remote studies of solar system comets (see Knight 
et al., in this volume). 

The similarity of 2I/Borisov to solar system comets ex- 
tended to the detection by Drahus et al. (ATEL 13549) of 
two outbursts on 2020 Mar. 4 and 8, brightening the comet 
by ~ 0.7 magnitudes. A sequence of HST imaging from 23 
March to 20 April to investigate possible nuclear fragmen- 
tation was reported by [Jewitt et al.|(2020b). These images 
showed a discrete nuclear secondary on 30 March, 26 days 
after the first outburst, that was not visible in images 4 days 
later. They posited this transient nature by the secondary 
consisting of one or more boulders ejected from the nucleus 
during the outburst, disrupting weeks later near the time of 
their observations due to rotational spin-up. 

The measured properties of 1I/“Oumuamua and 21/Borisov 
are summarized in Table[I] 


6.2.4. Composition and evolution 


The first characterizations of the coma of 2I/Borisov 
were obtained even before it was given a provisional des- 
ignation, with early optical photometry and spectra of the 
coma giving a red reflectance spectrum typical of normal 
comets 
and similar to 11/‘Oumuamua 
[et al.|2018). Gas emission via the bright CN(0-0) emission 
band at 388 nm was first detected on 2019 Sep. 20-24 by 


[Fitzsimmons et al.|(2019) and |de León et al.|(2019). Spec- 


troscopy and narrow-band photometry over the following 
month was reported by (2019). As with the 
earlier spectroscopy, they did not detect the C2(0-0) emis- 
sion band at 517 nm and concluded the comet was similar to 
carbon-chain depleted solar system comets. However weak 
detection of C2(0-0) was reported by several observers from 


Table 1: A summary of measured properties of interstellar objects 


Quantity 11/‘Oumuamua Ref? 21/Borisov Ref? 
Dynamical Properties 

Perihelion date Tp 2017 September 09.51 [1] 2019 December 08.55 [11] 
Perihelion distance q [au] 0.256 [1] 2.007 [11] 
Eccentricity e 1.201 [1] 3.356 [11] 
Earth close approach A [au] 0.162 [1] 1.937 [11] 
Incoming velocity? Və [km s71] 26.420 + 0.002 [2] 32.288 [11] 
Non-grav acceleration A; x 108 [au d7?] 27.90 [1] 7.09 [11] 
Non-grav acceleration A» x 10° [au d~?] 1.44 [1] -1.44 [11] 
Non-grav acceleration A3 x 10° [au d~?] 1.57 [1] 0.065 [11] 
Physical Properties 

Absolute magnitude Ay 22.4 + 0.04 [3] >19.1 [12] 
Albedo Pv 0.01 — 0.2 [4] -= 

Radius (for py=0.04) ryn [m] 110 [3] < 500 [12] 
Rotation state complex, long-axis mode [5] - 

Rotation period P [hr] 8.67 + 0.34 hr [5] —- 

Axis ratio a:b >6:1 [6] —- 

Spectral slope Sy [% per 100 nm] (7-23)+3 [6,7] 121 [8] 
H20 production Q(H20) [molec st] < 6.1x10?5 @ 0.38 au (obs) [9] 1.1 x 1027 @ 2.0 au (obs) [13] 
H20 production Q(H20) [molec s~'] 4.9 107° @ 1.4 au (model) [10] - 

Hyper volatile (CO?) Q(X) [molec s~!] 4.5 x102° @ 1.4 au (model) [10] - 

CO», production Q(CO2) [molec s~!] < 9 x10?? @ 2.0 au (obs) [4] - 

CO production Q(CO) [molec s~'] < 1 x10?4 @ 2.0 au (obs) [4] 4.4 x10? @ 2.0 au (obs) [14] 
Dust production Q(dust) [kg s7}] < 1078 @ 1.4 au (obs) [3,7] 2-— 50 @ 2.6 au (obs) [15,16] 


t vo is the velocity an object on a hyperbolic orbit has infinitely far from the central body (here, the Sun); ‘Reference Key: 


Jones et al.{2018); [3] (2017b); [4] the CO production rate is corrected from the number reported by 
018); [5] [Belton et al.| (2018); [6] (2019); [7] (2017); [8] [Jewitt and 
[Seligman] (2022); [9]|Hui and Knight\(2019); [10]|Micheli et al.|(2018); [11] JPL Horizons orbital solution #53, [12] Jewitt] 
a al (20200). 1131 King eral (2020, [14] Conner er al|(2020, T15] Jewin and Lan (9019), [161 de Leon eta OI), 
For production rates of other species for 21/Borisov, see the summary in|Jewitt and Seligman|(2022). 
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November (i.e. 2020), and by perihelion on 


8 December 2019 the C2/CN ratio was formally consis- 
tent with both depleted and normal alla chain cometary 
abundances (see Fig. 4 of eee ). These 
spectra and those presented by 1.) (2021) also 
showed other normal cometary molecular emission features 
due to species such as NH, NH2, and CH. Hence the optical 
species showed a significant evolution during this time as 
the heliocentric distance decreased from 2.7 au to 2.0 au 

The first constraint on the sublimation rates of primary 
ice species came from detection of the [OI] 6300 A line by 
(2020). Assuming this came solely from sub- 
limation and subsequent dissociation of water, they derived 
a production rate of Q(H2O) = (6.3 + 1.5) x 107° mol 
s~' on 29 October at r, = 2.38 au. (2020) 
used the UVOT telescope onboard the NASA Swift satel- 
lite to observe the OH (0-0) and (1-1) emission between 
280 — 320 nm, again resulting from dissociation of water 
molecules. The derived water production rates showed a 
peak of Q(H2O) = (1.1 + 0.1) x 1077 mol s~? near peri- 
helion, but with higher production rates pre-perihelion than 
post-perihelion. 

Given the significant activity at r, > 3 au, it was clear 
that other volatile species such as CO or CO% should be 


present in 21. (2020) used HST to measure 


the CO fluorescence bands at 140-170 nm and correspond- 
ing production rates during and after perihelion. 


(2020) used ALMA to measure CO (J = 3 — 2) 
emission along with HCN (J = 4 — 3). Both these 


teams showed that the CO abundance was exceptionally 
high in 2I/Borisov, being on par with the H2O abundance 
at Q(CO)/Q(H20) ~ 0.35 — 1.55. Only rare CO-rich 
comets such as C/2016 R2 had previously exhibited such 


CO-rich compositions (see (2019) and the 


chapter by Biver et al. in this volume). 

found that the CO production rate 
stayed relatively constant around perihelion, while that of 
H20 fell soon after, implying a significant change in the 
near-surface abundance ratios of these ices. This may have 
been due to seasonal effects of nucleus insolation, cou- 
pled with a non-heterogenous composition or erosion of 
the surface, similar to that seen by ESA Rosetta in some 
species at comet 67P/Churyumov-Gerasimenko 
jet_al.|[2020), although for that comet the CO/H20 ratio 
was relatively constant within 2.5 au. Similarly, as men- 
tioned above, the abundance ratios of the primary opti- 
cal species C2/CN also changed significantly as the comet 
moved towards perihelion. Finally, shortly after the dis- 
covery of significant neutral metal line emission in opti- 


= spectra of non-sungrazing comets (Hutsemékers et al. 
2021), these were also found in spectra of 2I/Borisov 


Guzik and Drahus\2021} |Opitom et al.J2021). The abun- 
Pare ratio of log(Q(Nil)/Q(Fel) 0.21 + 0.1 was 
within the range measured for solar system comets, as was 
Q(NiI + Fel)/Q(CO). 

Taken together, although several authors concluded that 
the optically active species were consistent with normal so- 
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lar system comet abundances near perihelion, these stud- 
ies point to a significantly evolving coma composition as 
2I/Borisov passed the Sun. It is unclear to what degree this 
evolution was caused by its interstellar nature, although the 
CO-rich composition is clearly unusual. 


6.3. 


With the discovery of 1I/‘Oumuamua and 21/Borisov, as- 
tronomers could finally estimate the true space density of 
ISOs in the Solar neighborhood, albeit with large uncer- 
tainties given N = 2. A complication also arises given 
their different apparent nature during their encounters, as 
to whether they may come from different populations of 
ISOs. Inert ISOs like 1/‘Oumuamua appear to have a lo- 
cal space density of ~ 0.1 to 0.3 au~3 at H = 22 (Meech 


2017b 2018). This approximately agrees 


an aa discovery upper limits for inert ISOs (Engelhardt 
et al.\2017) when scaled to this absolute magnitude. Some 


aoe reported a possible tension between the inferred 
mass density contained in ISOs and the galactic stellar den- 
sity, with 1J/‘Oumuamua implying an exceptional amount 
of mass-loss from planetary systems e.g. Horror 
However the 1) (2019) demon- 
strated that the space density of ~ 10! pe~? = 0.1 au7 
was compatible with possible size distributions. 

The space density of ISOs exhibiting normal cometary 
activity like 2I/Borisov is less well constrained, even though 
most earlier upper limits assumed cometary activity. |Engel-| 
[hardt et al.|(2017) found a pre-21/Borisov 90% upper limit 
of < 1.4 x 107 ~3 for a nuclear absolute magnitude 
H < 19, assuming ie entire sunward surface was active. 
used the statistics of LPC discoveries, 
together with assuming that ISOs share the same velocity 
distribution as stars within 100 pc, to derive a lower density 
of 7 x 1075 au~?. This is only a factor 2 less than {Engel-| 
[hardt et al.|(2017), and is plausible given the extra survey 
time since that study. 

If 11/‘Oumuamua and 21/Borisov come from the same 
ISO population, then Fig. [10|provides constraints on the lu- 
minosity distribution. We use absolute magnitudes to avoid 
any assumptions on albedo, although it is important to re- 
member these have been derived using assumed phase laws. 
If they share the same albedo then this also constrains the 
size distribution. For the interstellar space density from 
1]/‘Oumuamua we take the space density and confidence 


limits from|Levine et al.| (2021), and for 2I/Borisov we take 
the space density of |Eubanks et al.| (2021) together with 
Poisson uncertainty ranges calculated for a single detection 


at both magnitudes. 

We assume a single power-law relationship for the num- 
ber density against size, as with only 2 points a more com- 
plex relationship is unwarranted. An obvious size distri- 
bution to test is the theoretical differential size distribution 
for a population in collisional equilibrium, where the num- 
ber of objects at a radius R is given by dN/dR x R~4 with 
q = 3.5 (Dohnanyi\1969). A second size distribution would 
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Fig. 10.— Number densities for 1I/‘Oumuamua and 


2I/Borisov as a function of their absolute magnitudes. Also 
shown are differential brightness power-laws with expo- 
nents a = 0.5 for collisional equilibrium, and the range 
a = 0.46 to 0.2 measured for solar system comets. All are 
anchored to the absolute magnitude of 1I/‘Oumuamua. 


be that measured for solar system comets, normally given 
as a cumulative size distribution N(> R) œ R!~4. Several 


investigators have found q œ 2.9 (Meech et al.|2004 
[grass et al.|2011 2013). A more recent 


analysis of NEOWISE comet observations by [Bauer et al.| 
(2017) found a debiased value of q = 3.3 for JFCs similar 


to the earlier studies, but a shallower debiased g = 2.0 for 
LPCs. Note that in comparison with the mass distribution 
N(m) x m7? used in §4.2| the power-law exponents are 
related by q = (3p — 2). 

If we transform these size distributions to differential 
absolute magnitude distributions N(H) œ 10°”, the ex- 
ponents of the two forms are related by a = (q — 1)/5. 
The above size distributions imply a theoretical magnitude 
distribution with œa = 0.5, and measured values of a 
0.2 — 0.46. It is clear from Fig. [10] that 1J/‘Oumuamua 
and 2I/Borisov do not match these distributions. 

There are a number of possible explanations. First, they 
could be from different populations of ISOs, each of which 
has a more standard size distribution but with significantly 
different space densities. If they instead reflect the true 
size distribution of ISOs, this implies the size distribution 
is steeper than expected. As explained in ejection 
processes can produce an increase in the relative number 
of smaller bodies, consistent with Fig. An ejection- 
produced steepening of the size distribution would also ex- 
plain the mismatch with the observed LPC population that 
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originates from our Oort cloud, where objects have under- 
gone similar long-term exposure to the ISM but not ejec- 
tion. Finally, there could be a strong variation in albedo 
between bodies, possibly due to their individual histories 
in the ISM or, with 1I/‘Oumuamua, by its close perihelion 
passage pre-discovery. Identifying which of the these ex- 
planations is most likely will require further ISO detections 
for a much better measurement of the magnitude/size dis- 
tribution. 


6.4. CNEOS 2014-01-08 


While 1I/‘Oumuamua and 2I/Borisov are macroscopic 
ISOs, the existence of interstellar dust particles — micro- 
scopic ISOs — has been known for some time from space- 
craft detection (see the review by|Sterken et al.|2019). The 
Advanced Meteor Orbit Radar multi-station complex has 
provided extensive data on interplanetary dust characteris- 
tics. The dynamical information that is collected allows 
identification of discrete source regions. One of the main 
sources appears to be the 8 Pic disk (Baggaley|2000). Dur- 
ing its cruise phase, the Stardust mission captured particles 


from the oncoming interstellar dust stream (Westphal et al. 
2014). 


In between dust and ISO sizes should lie objects ~1 cm 
-10 m in size. Fireballs produced by objects at these sizes 
are regularly detected via US DoD satellites, which mea- 
sure both the luminosity of the fireball and (above a lumi- 
nosity limit) the entry velocity vector. A search through 
the publicly accessible fireball catalogue led|Siraj and Loeb| 
to identify a potential interstellar impactor, CNEOS 
2014-01-08. This object entered the Earths’ atmosphere 
over the Western Pacific with a reported entry velocity of 
44.8 km/s. Although formal measurement uncertainties are 
not published, DoD personnel communicated that “the ve- 
locity estimate reported to NASA is sufficiently accurate to 
indicate an interstellar trajectory”. When integrated back- 
wards, the reported velocity vector implied an original orbit 
with e = 2.4 + 0.3, hyperbolic at the 3-ø level. 

While the identification of an interstellar origin appears 
sound based on the available satellite data, there is a ten- 
sion with the fireball and meteor data regularly obtained 
by ground-based optical and radar meteor surveys. 
identified 2 hyperbolic orbits out of 1739 me- 
teors observed with the Canadian Automated Meteor Ob- 
servatory, but careful inspection ruled them out as measure- 
ment error. Similarly, in an analysis of 824 fireballs de- 
tected by the European Fireball Network, 
identified two objects that had orbital eccentricities 
e > 1 at the 3-o level, but the absolute values were near 
unity and they concluded there was no strong evidence for 
a hyperbolic nature. In an analysis of nearly 4000 meteors 
observed by the FRIPON project, re- 
port an upper limit of interstellar meteors of 0.1% but sus- 
pect the true value to be much lower. Collating several me- 
teor orbit catalogues derived from photographic, video and 


radar systems,|Hajdukova et al.|(2020) concluded that there 


was no convincing evidence of any interstellar meteors be- 
ing detected in 25 years of data covering many thousands 
of objects. That said, a careful analysis of 4 x 10° meteors 
detected by the Canadian Meteor Orbit Radar led {Froncisz| 
to identify 5 possible interstellar meteors. We 
note that this possible interstellar fraction of ~1 in 80,000 
meteors is significantly lower than the ~1 in 900 implied 
for the current CNEOS fireball database. 

It is clear from these latter studies that measurement un- 
certainties are extremely important in interstellar meteor 
identification, due to the short timespan over which the ob- 
servations are obtained. Until such a time that a quantitative 
description of the uncertainties on the DoD satellite mea- 
surements is forthcoming (and indeed their actual values), 
some uncertainty unfortunately still remains concerning an 
interstellar origin for CNEOS 2014-01-08. It therefore fol- 
lows that the pursuit of interstellar meteor detections with 
quantifiable uncertainties is a highly worthwhile endeavor. 


6.5. 


As reported in the Planetary Decadal survey 


Academies of Sciences and Medicine\|2022), as of 2021 
about 95% of the NEOs greater than | km in diameter have 


been found. These objects are capable of causing global 
effects upon impact. Objects that are larger than 140 m in 
diameter can cause regional destruction and to date, only 
about 33% of these have been discovered. Estimates show 
that there might be more than 10° objects >50 m in diame- 
ter which could cause destruction on urban scales, and less 
than 2% of these have been found. In 2005 Congress di- 
rected NASA to find 90% of all NEOs larger than 140m, 
and to date this directive has not been met. Once it is met, 
the threat from impacts from near earth asteroids will be 
minimized, but this will not address the risk from objects 
on long period Oort cloud comet (OCC) trajectories. It is 
likely that OCCs are responsible for most of the very large 


impacts on Earth (Jeffers et al.|2001). 
The typical NEO encounter velocity peaks around 15 


km/s, with a range up to 40 km/s (Heinze et al.|/2021), 
however the encounter velocities for OCCs will typically 
be around 54 km/s ranging from 16 up to 72 km/s 
[et al.2001). ISO trajectories will be very similar to those 
of OCCs, but the encounter velocities can be even higher. 
At the time of close approach to the Earth on 2017 Oct. 
14, the relative velocity of 1]/‘Oumuamua with respect to 
the Earth was 60.3 km/s. It passed within 63 Earth-moon 
distances and we did not and could not have discovered 
it until after it had passed. Objects like 1I/‘Oumuamua 
which are on OCC-like orbits but which do not have any 
detectable activity, have small nuclei, and low albedo (such 


as the Manx comets, see|Meech et al.|2016) are particularly 
dangerous because they will be harder to detect. 


also noted that there is a larger bias than ex- 
pected against finding small fast moving objects because of 
the streak length on the detector, so inactive ISOs will be 
particularly hard to find and we are more likely to be miss- 
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ing more of these. 

These OCC and ISO trajectories will be distributed 
across a wide range of inclinations. As discussed in 
the powerful new Vera C. Rubin Observatory (Rubin) 
Legacy Survey of Space and Time (LSST) will begin in 
mid-2025. There is a concerted effort to optimize the sur- 
vey effort to benefit solar system science, including de- 
tection of OCCs, NEOs and potentially hazardous objects 
(Schwamb et al.2023). 11/‘Oumuamua was discovered at 
1.22 au moving at a rate of 6.2° per day, i.e. typical of a 
faint nearby NEO that the LSST will be optimized to detect. 


7. THE NEXT DECADE OF EXOCOMET AND ISO 
STUDIES - WHAT DON’T WE KNOW? 


7.1. 
7.1.1. 


Exocomet Systems 
Probing outer exocomet populations 


Current and under-construction facilities should lead to 
significant advances in the field of exocomet science. An 
important area of study will be to understand the dynami- 
cal evolution of exocomets in planetary systems, which re- 
quires both high angular resolution and sensitivity at dust- 
emission wavelengths. Near to mid-IR dust imaging will 
take place with the 6.5m James Webb Space Telescope 
(JWST) , the 39-m ESO Extremely Large Telescope (ELT) 


with its MICADO and METIS cameras (Brand et al.\2010 


[Davies et al.2018), and large, deep ALMA surveys such 
as the ongoing ARKS (The ALMA survey to Resolve ex- 


oKuiper belt Substructures, Marino et al. in prep.). These 
will deliver unprecedentedly detailed images of dust emis- 
sion from exocometary belts, revealing planet-belt inter- 
action and allowing us to infer the dynamical fate of ex- 
ocomets within belts, be it ejection, inward scattering, or 
survival within the belt. 

Gas observations offer a unique window into exo- 
cometary compositions in young planetary systems just 
after formation. Near- to mid-IR spectroscopy with JWST 
combined with ground-based near-IR high-resolution spec- 
trographs (like the newly-installed CRIRES+ instrumen- 
tation on the ESO Very Large Telescope and later the 
ELT+HARMOND will be used to probe ro-vibrational tran- 
sitions of new molecules undetectable at mm wavelengths, 
such as CO2, CH; and of course H2O and its photodis- 
sociation product OH. This will constrain exocometary 
molecular gas compositions within outer belts, with the 
potential to confirm the gas origin, determine the physical 
release mechanism, and ultimately the ice composition of 
young exocomets - a missing evolutionary link between 
outer protoplanetary disks and solar system cometary/KBO 
compositions. 

In addition, ALMA CO and C I line imaging will con- 
strain the radial and vertical distribution and hence the evo- 
lution of gas within exocometary belts; this will be com- 
bined with detail kinematic analysis, thanks to ALMA’s 
maximum spectral resolution, to reveal the environmental 
conditions (temperature, bulk density) of the gas, and eval- 


uate the dynamical influence of planets on the evolution of 
exocometary gas. 


7.1.2. Probing the inward transport of exocomets 


The highest priority for the future of exocomets in the 
inner regions should be to expand the number of detections 
and detected systems - currently dominated by 8 Pic - and 
study their compositions and dynamics as they enter the in- 
nermost regions of a planetary system. This will enable us 
to probe their formation location and the mechanism lead- 
ing to their inward transport - with the potential of linking 
their origin to detected exoplanets, and imaged populations 
of exocomets in outer belts. 

JWST is especially sensitive to warm dust emission in 
the interior regions of nearby systems, either through its di- 
rect imaging or aperture masking modes, and this can be 
used to set limits on inward scattering and transport of icy 
exocomets from the outer belt towards the terrestrial re- 
gion (Marino et al.|2018a). The Roman Space Telescope 
is a 2.4-m wide-field optical and near-IR facility, currently 
planned to launch on 2027. Its coronagraph instrument is 
also expected to expand the number of systems with exozo- 
diacal light detections in the terrestrial planet region 
glas et al./2022). Imaging of exozodiacal dust may for the 
first time also be possible around the nearest stars with the 
ELT (Roberge et al.|2012). 

The identification and confirmation of more 8 Pic-like 
systems requires high resolution echelle spectrographs with 
large aperture telescopes to provide high signal-to-noise 
per wavelength bin. Observations with new high-resolution 


spectrographs on the ELT such as ANDES 
[et _al.|[2013) and METIS (Brandl et al.|2010) will exploit 


the ELT’s large collecting area with their high spectral 
(R~100,000) resolution and wide wavelength coverage. 
This will allow more sensitive searches for absorption in 
edge-on systems, both from narrow, stable gas at the stellar 
velocity and variable red- and blue-shifted features from ex- 
ocometary gas. Monitoring observations looking for high- 
velocity exocomet features should aim to cover volatiles as 
well as metallic atomic species across a broad wavelength 
range, as one of the highest priorities is to understand the 
nature and volatile content of star-grazing exocomets. 

Finally, continued exploration of TESS and CHEOPS 
stellar light-curves to look for asymmetric exocomet tran- 
sits are likely to lead to new discoveries. These will be fur- 
ther supported by the PLATO transiting exoplanet mission 
(Rauer et al.|2014), with its factor ~ 5 increased photo- 
metric accuracy over TESS, currently scheduled for launch 
in 2026. This will contribute to understanding the occur- 
rence rate of star-grazing exocomets around large samples 
of stars, pushing towards smaller transits. 


7.2. Galactic evolution of ISOs 
7.2.1. Galactic effects on ISOs 


In $4.2] we saw that ISOs experience various influences 
from the environment in the form of thermal heating and ex- 
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posure to radiation and particles of various kinds. The ques- 
tion is whether ISOs also affect the environments they are 
passing through. The discovery of II/‘Oumuamua taught 
us that the ISM contains not only gas and dust but that ISOs 
are a natural third component of the ISM. 

ISOs also pass through molecular clouds (see section 
4.3), with young ISOs more likely to make such an expe- 
rience than older ones (Pfalzner et al.J2020). Such a cloud 
passage might alter not only the ISO surface and structure 
but also the direction of their path, while at the same time 
decelerating or accelerating them. Depending on the size of 
the actual molecular cloud 101° — 107° ISOs reside in them 
at any time. Being so small in mass, their individual effect 
is negligible, but it could be different if co-moving streams 
of ISOs pass through molecular clouds. 

However, although every ISO travels through the ISM, 
certainly not every ISO passes through a molecular cloud 
and even a protoplanetary disk. However, the probability of 
ISOs passing through molecular clouds is surprisingly high. 
In the solar neighborhood, an ISO spends ~ 0.1%-0.2% 


of its journey passing through molecular clouds (Pfalzner| 
[et _al.|/2020). This value increases for young ISOs close 


to the Galactic center. In comparison, passing through an 
existing protoplanetary disk is a much rarer event due to 
its much smaller cross-section. Nevertheless, {Grishin et al.| 
estimate that at least 104 ISOs larger than 1 km cross 
any protoplanetary disk around field stars, with the number 
increasing to 10° for star cluster environments. 

Under certain conditions, molecular clouds collapse and 
form stars. It seems that ISOs take part in this process 
to some extent. This would mean that ISO could become 
concentrated in collapsing molecular clouds 
[2021b). However, these simulations indicate that there is a 
competitive capture process at work that favors the capture 
of ISOs by massive star clusters. 

In the current planet formation scenario, there exist some 
difficulties in proceeding from m-sized to km-sized objects 
(see (2.3). It has been suggested that ISOs might easily 
overcome the meter-sized growth barrier by acting as seeds 
to catalyze planet formation. Two different scenarios have 
been suggested. theorize that ISOs 
could be captured from the local star formation region into 
the disks surrounding young stars . Alternatively, [Pfalzner| 
propose the ISOs in the ISM taking 


part in molecular cloud collapse would become a natural 
component of the forming disk, without the need for addi- 
tional capture. They give a first conservative estimate of the 
order 10'! ISOs typically being incorporated into forming 
star-disk systems. (2022) found 
a similar total number when considering metre-scale ISOs 
and larger, but also noted the number in the disk could in- 
crease by 2-3 orders of magnitude in cluster environments. 
If ISOs take part in cloud collapse leading to star forma- 
tion, this would also mean that ISOs are incorporated into 
forming stars. The question is how many ISOs end up in 
stars? 

When ISOs pass through molecular clouds and are cap- 
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tured in planet-forming disks, they can be affected by these 
environments. Possible mechanisms are devolatilization 
and erosion. Both will lead to changes in the ISO size dis- 
tribution. Whether the ISOs are altered depends primarily 
on the gas/dust density, which increases from an ISM envi- 
ronment to molecular clouds and protoplanetary disks. 

Finally, (2018) point out that exo-Oort clouds 
around stars that produce supernovae will be irradiated and 
lose surface volatiles, but the exocomets there will survive. 
As they then drift away from the supernova remnant they 
will form a natural population of devolatilized ISOs, pos- 
sibly somewhat similar to 1I/‘Oumuamua. However super- 
novae are rare events, and such ISOs should not form a ma- 
jor component of the population. 


7.2.2. ISO effects on the galaxy 


Over the last four years, it has become apparent that 
ISOs are a natural component of molecular clouds. The 
question is whether the presence of the huge number of 
low-mass ISOs does “contaminate” star forming regions. 
This depends on the mass of the ISOs present in molecu- 
lar clouds, but also on the different chemistry between star 
forming systems and on whether ISOs are evenly spread or 
come in streams from their parent systems. 

So far this question has been studied more from the per- 
spective of individual already forming planetary systems. 
This question of whether material can be spread from one 
planetary system to another is of long standing interest (Val- 
tonen and Innanen] 1982, 
Brasser et al./2006). More recently, estimates of the 
likelihood of ISO capture events have been performed for 
the early phase in young star-forming clusters and the local 


galactic neighborhood (Hands et al.\2019 


2021}|Napier et al.|2021). In young clusters the stellar den- 
sity is much higher than in the field, therefore, the ISO 


capture rate per Myr is naturally much higher. However, 
whether the capture probability is in general higher during 
the short cluster phase than over the entire lifetime of a star 
depends strongly on the assumed cluster properties and the 
stellar mass. The latter determines the capture cross sec- 
tion and the lifetime of a star. Besides, GAIA data strongly 
indicate that clusters seem not completely dissolved within 
10 Myr, but move together in an unbound state for a con- 
siderable length of time. Thus assuming field densities for 
the stars and the ISOs underestimate the capture rates in the 
first few hundred Myr. This requires further investigation in 
the future. 

On average every star will contribute ISOs to the galactic 
population, and this material might be injected into the in- 
ner regions of star systems. suggest 
that this material can impact both the host stars and their 
planets, enriching their atmospheres, and that understand- 
ing the composition of the ISO population will have impli- 
cations for post-formation exoplanet atmospheric composi- 
tion. 

Currently there exists no definitive evidence that any 
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gravitationally bound objects in the solar system are of ex- 
trinsic origin (Morbidelli et al.[2020). However, the discov- 
ery of 11/‘Oumuamua and 2I/Borisov sparked wide-ranging 
speculation regarding the possibility that our solar system is 
being more broadly contaminated by minor bodies of extra- 


solar origin (Siraj and Loeb\|2019 
2020}|Hands and Dehnen\2020). There exists a wide range 


in the estimated masses of planetesimals accreted from 
other stars while the Sun lived within its birth cluster from 


1075 Mg to a third of the Oort cloud population. 


et al.|(2021) finds that about 10° more ISOs were captured 
p 


during the cluster phase than accreted subsequently from 
the field, while the current steady state factor is 10* and the 
total mass of surviving captured ISOs is ~ 107° Mg. To 
put this into perspective, the mass of the Oort cloud is es- 


timated to range from 0.75 + 0.25 Mg (Brasser|2008) to 
1.9 Me (Weissman|1996). With a typical comet-mass of a 
few times 10°* to 10°* kg (Sosa and Ferndndez}2011) the 


Oort cloud contains approximately 10°°-10°* comet-sized 
objects. Thus the captured ISO population is probably ex- 
tremely small in comparison to the primordial comet popu- 
lation in the Oort cloud. 

Current capture of ISOs into the solar system has also 
been explored. suggest that plan- 
ets capture ~2 ISOs every 1000 years which would re- 
sult in 8 ISOs captured within 5 au of the Sun at any 
time. A perennial candidate for a captured ISO is comet 
96P/Machholz, due to its highly anomalous optical spec- 
trum that is dominated by NH2, with strong depletions of 


carbon-based molecules (Langland-Shula and Smith\2007 
Schleicher|\2008). Its orbital elements of q 0.124, 


i = 58° and e = 0.959 mark it as one of the lowest per- 
ihelia comets known, although its high inclination is due to 
strong Kozai-Lidov oscillations (McIntosh||1990). Its nu- 
cleus was studied by Sr oR compare with 
11/‘Oumuamua to search for possible effects of small peri- 
helion passages. The strongest evidence for it being a cap- 
tured ISO remains its peculiar coma composition, but this 
is not a definitive marker. 


7.2.3. 


The existence of ISOs in the Milky Way makes it likely 
that other galaxies are also populated with large numbers of 
ISOs. One can also speculate about the number of objects 
that have been ejected from their parent galaxy and become 
Intergalactic Objects. These would clearly be difficult to 
create, given the huge velocities required; at the Sun’s posi- 
tion the escape velocity from the Milky Way is ~ 500 km/s 


(Koppelman and Helmi|2021). While some tens of hyper- 


velocity stars escaping the Milky Way have now been iden- 
tified (Li et al.[2021), most are thought to have attained their 
velocities via encounters with massive black holes, whose 
dynamical and radiation environments are hardly conducive 
for ISO survival. 

However, we note that while escape to intergalactic 
space is highly unlikely for an ISO, we can let the galaxy 


Intergalactic Objects 


carry it to a new home. Galactic mergers have taken place 
throughout the history of the Universe, with 5 minor merg- 
ers identified for the Milky Way (Kruijssen et al.|2020). An 
unknown fraction of the ISOs in those galaxies would have 
merged with the galactic population, and of course sub- 
sumed stellar systems would have ejected their own ISOs 
through the usual methods described in Hence given 
enough time as ISO detections increase, there is a finite but 
as-yet unquantified probability of detecting an ISO carried 
here by its own galaxy. 


7.3. The range of ISO properties 


7.3.1. Future detections and surveys 


The trajectory of 1I/‘Oumuamua brought it in from 
above the plane of the solar system from the direction of 
the constellation Lyra (Fig. m. It came to perihelion in- 
side the orbit of Mercury at q = 0.255 au on 2017 Sep. 9. 
and was discovered post-perihelion at r = 1.22 au after it 
had made its close approach to Earth on 2017 Oct. 14. This 
close approach was at 0.162 au, passing within 63 lunar dis- 
tances. Images from the Catalina Sky Survey detected it on 
Oct. 14, but it was noted only after it had been discovered. 


| 1837 — Passed inside 1000 au 
Jan 18, 2017 — inside 5.2 au 
Aug 10, 2017 — inside 1.0 au 


May 3, 2018 — outside 5.2 au 
Jun 2022 -30 au 
Feb 2024-39 au 
Dec 2025-50 au 
2038-121 au 
2196 — 1000 au 


Sep 9, 2017 - perihelion q = 0.255 au 
Oct 11, 2017 - outside 1.0 au 
Oct 14, 2017 — Earth approach A = 0.162 ay 


18 Oct — 2 Nov 2017 


A/2017 Ud path 


Comet orbit * 


Fig. 11.— Path of 1]/“Oumuamua as it entered the solar 
system illustrating why it was not discovered sooner. 


11/“Oumuamua passed through the SOHO and STEREO 
fields of view in early 2017 Sep., near perihelion 
2019). Because of the large phase angle and the 


extreme forward scattering, this would have enhanced the 
brightness of any dust around the nucleus by ~ two or- 
ders of magnitude but neither satellite detected it. It could 
not have been discovered any earlier than it was because 
11/“Oumuamua was just coming out of solar conjunction 
as seen from the ground and was faint (mag~23.7 at 0.79 
au on Oct. 2). Prior to moving into solar conjunction in 
early August it would have been fainter than mag~25 - both 
times too faint to be discovered by any of the existing sur- 
veys. 1I/‘Oumuamua was inside the orbit of Jupiter for less 
than 1.3 yrs. The second ISO was discovered by an amateur 
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searching close to the horizon in twilight — something that 
the large telescope surveys cannot do (see Bauer et al., this 
volume). 

The Rubin observatory in Chile has an 8.3m survey tele- 
scope with a camera that covers a FOV of 9.6 deg?. It 
will begin regular survey operations in mid-2025 that will 
provide a much greater depth than all previous surveys, 
reaching mag~24.7 and scanning the accessible sky every 
3 days (Jones er al 2009) (2022) have es- 
timated the number of discoveries that might be made by 
the LSST by generating a synthetic population of ISOs, as- 
suming no cometary activity. They predict that the survey 
will find on the order of 1 ISO per year, but the number 
can be larger than 1 when cometary activity is considered. 
Ironically, however, even had the LSST survey been op- 
erational in 2017, it would have been unlikely to discover 
1I1/ʻOumuamua because the ISO would have had an aver- 
age brightness brighter than mag~24.7 for less than a week 
before moving into solar conjunction in early 2017 Aug. 


7.3.2. Probing origins and evolution 


The standard mechanisms of ISO creation via ejection of 
exocomets from their home systems are described in detail 
in 1I1/ʻOumuamua’s discovery resulted in the proposal 
of additional possible creation routes. Many of the forma- 
tion mechanisms have been suggested in an attempt to ex- 
plain the origin of the unusual shape for 1I/‘Oumuamua. 
These include fluidization to a Jacobi ellipsoid during the 
red giant phase of a star 2018), interstellar ablation 


(Vavilov and Medvedev\2019), collisional elongation (Sug- 
iura et al.||2019), formation as an icy fractal aggregate 


2019a) and a tidally disrupted planetesimal 
that passed close to growing giant planets, from which the 


volatiles were removed during a close stellar passage 
Zhang and Lin}2020). 

Two 1I/‘Oumuamua origin scenarios suggested a volatile 
composition based on homonuclear diatomic molecules 
with no dipole moments and no vibrational spectral lines, 
which would explain why no gas was detected. In an 
attempt to both explain the shape and acceleration of 
11/‘Oumuamua, suggested 
that it was composed of molecular hydrogen ice. In their 
model, mass wasting from sublimation far out in the solar 
system could explain both the shape and the acceleration. 
They proposed that the ISO formed in a cold dense molecu- 
lar cloud core. The freezing point of Hə is around 14K, and 
the lowest temperature cloud cores are around 10K. Hə has 
no dipole moment and can only be detected in the far-UV 
or in the infrared through its rotational lines from space. 
An alternate suggestion was made by [Desch and Jackson] 
who proposed that 1I/‘Oumuamua was a collisional 
fragment from an exo-Kuiper belt. Many large KBOs in our 
solar system exhibit Nə ice on their surfaces. Gaseous Nə 
has no bending mode vibration spectral lines and is infra- 
red inactive. Its spectral features are in the UV where no 
observations were taken. This makes it an attractive volatile 


to explain how 1I/‘Oumuamua could have unobserved out- 
gassing. suggests that neither scenario 
is plausible because large H2-ice bodies are not likely to 
form in cloud cores and that the impacts in a Kuiper belt 
are not sufficient to generate large No fragments. 

In order to reconcile some of the observational and 
model inconsistencies, pro- 
pose that the acceleration of 1I/“Oumuamua is caused by 
release of molecular hydrogen that formed by cosmic ray 
processing of water ice. These authors argue that labora- 
tory experiments show that Hə is efficiently produced in 
water ice during radiation by high energy particles. They 
propose that the Hə gas is released over a wide range of 
temperatures during annealing of amorphous water ice and 
that there is sufficient gas released to account for the ob- 
served non-gravitational acceleration. 

Finally, regarding origin location within a disk, |Seligman| 
propose that measuring the C/O ratio of an 
ISO can be a tracer of whether the ISO formed inside or 
outside the ice line in its home star system. 


7.4. 


1I‘/Oumuamua was accessible to ground based tele- 
scopes for less than a month, a little longer using space fa- 
cilities. After this brief period of observation it was found 
that the characteristics were quite different from what was 
expected from the first ISO, and this one discovery has en- 
ergized a new interdisciplinary awareness in the study of 
planet formation. However, a more detailed investigation 
of ISOs with an in-situ mission presents unique challenges: 
the orbits may have high inclinations and the encounter 
speeds are typically high (10s of km/s) which means a short 
encounter. There are increased risks with high velocity fly- 
bys if the ISO ejects large dust particles. The Giotto space- 
craft was destabilized with a cm-sized impactor at a relative 
velocity of 68 km/s (Bird et al.J1988). 

NASA’s competed mission calls are not compatible with 
missions that are responsive to new discoveries, i.e., mis- 
sions that do not have a target at the time they are pro- 
posed. This is relevant for ISOs as well as OCCs, Manx 
comets, and hazardous NEOs. Two approaches have been 
suggested to explore these targets: spacecraft in storage, 
ready to launch following target discovery and spacecraft in 
a standby orbit, as is being done by ESA’s Comet Intercep- 


tor mission (see Snodgrass et al this volume, and 
2019 2021). Typically, the target 


will be unknown at the time of spacecraft development and 
this has an effect on the definition of basic spacecraft capa- 
bility (e.g., Av) and on the payload. Launch following the 
discovery of an ISO offers a greater flexibility in terms of 
target access but requires a very fast turnaround of a launch 
vehicle. A spacecraft in a standby orbit is more responsive 
to a target but has a more limited target accessibility. 

Using known OCC orbits, an assessment of the phase 
space of targets with available launch energies, maximum 
encounter speed and for times of flight of less than 10 years, 
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launching while the comet is more than 0.5 years from com- 
ing to perihelion, shows that for low launch energies, C3, 
there are almost no targets. (In astrodynamics, the charac- 
teristic energy C3 km?/s? is the measure of the excess spe- 
cific energy over that required to just barely escape a mas- 
sive body; C3 = v2.) There are also almost no targets unless 
the encounter speeds are greater than 10 km/s. This sug- 
gests that at present only fast flyby missions are an option 
for ISOs. For a comprehensive mission, this likely requires 
significant advancements in payload capabilities (e.g. small 
deployable satellites and autonomous navigation) 
[2022). 


There were several concepts developed for missions that 
could reach 1I/‘Oumuamua. (2018) 
estimated that for launch-on-detection scenarios, the wait 
time would be on the order of 10 years for a favorable mis- 
sion opportunity. (2020) explored a range 
of flyby trajectories which could launch in the early 2030’s 
and deliver a spacecraft to 1I/‘Oumuamua with relative ve- 
locities between 15-20 km/s, arriving between 2048-2052. 
Finally, {Miller et al.|(2022) proposed high-performance so- 
lar sail scenarios that would allow for a rendezvous with an 
ISO by maintaining a high potential energy position until 
the detection of an ISO and then matching velocity through 
a controlled fall toward the Sun. 

The LSST survey will increase the number of discov- 
eries by going much fainter, and therefore to larger dis- 
covery distances. As shown by (2017), 
ISOs are more common at larger perihelion distances. The 
fainter limiting survey magnitude would enable detection 
of a 1 km inactive nucleus with 4% albedo out to 5.5 
au, active ISOs much further. Although most detected 
ISOs may be too distant to easily reach with spacecraft, 
11/‘Oumuamua and 21/Borisov demonstrate ISOs exist that 
come within 2 au, and detection of an incoming ISO at large 
distances may provide enough warning time for either the 
storage or standby concepts discussed above. 


8. SUMMARY 


From the multitude of studies described in this chapter, it 
is clear that exocomets should be and are common in nearby 
stellar systems. Together with the physical mechanisms by 
which they can be lost to interstellar space, the local stellar 
neighborhood should be rich in ISOs, a prediction finally 
confirmed through the discovery of 1I/‘Oumuamua and 
2I/Borisov. So the question remains, how is 1I/‘Oumuamua 
not visibly active while 2I/Borisov is strongly active; is this 
due to origin, age, or the unseen 1I/‘Oumuamua perihelion 
passage? From models of our solar system (Shannon et al. 
(2015), some fraction of ISOs will not be cometary but more 
like C- or S-type asteroids. Based on the observation of 
Manx comets (Meech et al.\|2016), we should be able to 
distinguish asteroidal objects from inert comets. Is there 
a continuum of properties for ISOs, or have we already 
found two separate populations arriving from different ori- 
gin mechanisms? This sample of two very different ISOs 
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12.— Genetic relationships of early solar system planetesimals both in our solar system and from exoplanetary 


systems. We are just at the beginning of exploring the ISO population. Objects on OCC orbits, including ISOs, represent 
the largest reservoir of objects that has not been explored with an in-situ mission. HTC = Halley Type comet, JFC = Jupiter 


Family Comet). 


makes it difficult to predict what will 3I be like - will it 
be like planetesimals from our solar system or something 
different (see Fig. [12)? Given the increasing sensitivity 
of sky surveys, the current and next generation of opti- 
cal, IR and sub-mm facilities, plus the theoretical advances 
in cometary/exocomet structure and evolution, these open 
questions have a chance of being answered in Comets IV. 
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